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ABSTRACT 


Work  is  continuing  primarily  in  gas  phase  turbulent  mixing  and 
chemical  reactions  with  important  ' j  extensions  to  compressible 
(supersonic)  shear  layers..  This  work  is  divided  into  an  experimental 
program,  a  theoretical  and  modelling  program,  a  computational  program 
and  a  diagnostics,  instrumentation  and  data  acquisition  program.  In  the 
gas  phase  shear  layer  work,  ourv  investigations/  are  concentrating  on 
subsonic  shear  layer  free  stream  density  ratio  effects,  and  a  design 
effort  in  support  of  the  planned  extension  of  the  hydrogen-fluorine 
shear  flow  facility  to  supersonic  flows.  The  work  on  chemically 
reacting  liquid  shear  layers  has  been  published  in  its  final  form 
(Koochesfahani  &  Dimotakis  1986).  In  jet  flows,  ihe  measurements  of  gas 
phase  jet  mixing,  using  laser  Rayleigh  scattering  techniques  developed 
for  conserved  scalar  measurements  down  to  diffusion  space  and  time 
scales,  are  in  progress.  A  first  publication  has  just  appeared  in 
Combustion  &  Flame  (Miake-Lye  &  Toner  1987)  .on  an  experiment  ,  completed 
with  joint  support  from  the  Gas  Research  Institute,  -in  which  digital 
imaging  of  soot  in  turbulent  flames  was  used  to  describe  cthe  combustion 
flame  sheets  in  methane  flames. 

Analytical/computational  modeling  efforts  during  this  last  year 
have  included  the  *  development  of  a  quantitative  description  of  turbulent 
jet  mixing  and  chemical  reactions  including  finite  •  Damkohler*  number 
effects,’  supersonic  shear  layer  combustion  studies  of  finite  kinetic 
rate  (Damkphler  number)  effects  for  the  H2/F2/NO  and  i^/air  systems,*  a 
new  analytical  model  for  turbulent  shear  layer  mixing  and  chemical 
reactions,' , extensions  of  hydrodynamic  stability  calculations  to  include 
Mach  number  effects  in  supersonic  shear  layers.  Finally,  computational 
fluid  dynamics  studies  have  also  proceeded  using  vortex  methods  as  well 
as  more  fundamental  studies  using  the  machinery  of  non-linear  systems 
dynamics. 
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1.0  INTRODUCTION 


This  project  is  concerned  with  mixing  and  chemical  reactions  in 
moderate  to  high  Reynolds  number,  turbulent  shear  flows  whose  scope  has 
been  extended  during  this  last  year  to  cover  the  effects  of 
compressibility  on  turbulence  and  turbulent  mixing.  It  is  directed  and 
making  significant  progress  towards  addressing  fundamental  deficiencies 
in  the  conventional  treatment  and  understanding  of  this  subject. 

The  program  is  comprised  of  several  complementary  parts.  In 

particular , 

1.  an  experimental  effort, 

2.  an  analytical/computational  modeling  effort, 

3.  a  computational  fluid  dynamics  effort, 


and 


4.  a  diagnostics/data-acquisition  effort,  the  latter  as  dictated 
by  specific  needs  of  our  experimental  program. 

Our  approach  is  to  carry  out  a  series  of  detailed  studies  in  two 
well  defined,  fundamentally  important  flow  fields:  free  shear  layers 
and  axisymmetric  jets.  The  philosophy  in  effect  being  that  these  two 
flows  encompass  an  important  part  of  the  generic  properties  and  behavior 
of  these  flows.  If  the  behavior  in  more  complicated  geometries  is  ever 
to  be  mastered,  understanding  of  these  two  canonical  flows  needs  to  be 
recognized  as  an  important  milestone  and  proving  ground  of  proposed 
models. 

To  elucidate  molecular  transport  effects,  experiments  and  theory 
concern  themselves  with  both  liquids  and  gases.  A  review  of  these 
issues  appeared  in  the  AIAA  J.  last  June  (Broadwell  A  Dimotakis  1986), 
which  received  a  citation  from  the  editorial  staff  of  the  AIAA  J. 


3 


1 -Jun-87 


AFOSR-83-021 3  1987  Annual  Report 


The  gas  phase  studies  are  presently  being  extended  to  cover  Mach 
number  effects  in  shear  layer  mixing,  i.e.  the  effects  of 
compressibility  in  supersonic,  high  Reynolds  number,  reacting  and 
non-reacting  shear  layers.  Progress  in  these  areas,  realized  since  our 
last  annual  report  (Dimotakis,  Broadwell  &  Leonard  1986),  is  outlined 
below. 


2.0  SHEAR  LAYER  MIXING  &  CHEMICAL  REACTIONS 


Turbulent,  two-dimensional  shear  layer  mixing  and  combustion  is  of 
primary  interest  in  this  effort.  Studies  in  this  area,  published  under 
the  sponsorship  of  this  program  prior  to  the  beginning  of  this  reporting 
period  (15  April  1986  -  15  April  1987),  have  concerned  themselves  with: 


1.  a  simple  model  for  turbulent  shear  flow  mixing  (Broadwell  & 
Breideenthal  1982), 


2.  gas-phase  mixing  at  low  heat  release  (Mungal  1983,  Mungal, 
Dimotakis  &  Broadwell  1984,  Mungal  &  Dimotakis  1984),  including 
Reynolds  number  effects  (Mungal,  Hermanson  &  Dimotakis  1985), 
and  Damkohler  number  effects  (Mungal  &  Frieler  1985), 


3.  a  simple  model  for  the  entrainment  ratio  in  a  subsonic, 
two-dimensional  shear  layer  (Dimotakis  1986), 


4.  gas-phase  mixing  at  high  heat  release,  including  the  effects  of 
favorable  pressure  gradients  (Hermanson  1985), 


and 


5.  liquid-phase  mixing  and  chemical  reactions  (Koochesfahani  1984, 
Koochesfahani ,  Dimotakis  &  Broadwell  1985,  Koochesfahani  & 
Dimotakis  1985) 


Work  published  or  in  progress  during  the  current  reporting  period  is 
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outlined  below. 


2.1  Heat  release  effects 

A  paper  (Hermanson,  Mungal  &  Dimotakis  1987)  has  just  appeared 
( AIAA  J. ) ,  documenting  part  of  the  progress  in  this  area.  A  first  draft 
of  the  manuscript  for  the  major  publication  derived  from  this  effort 
(Hermanson  1985).  intended  for  submission  to  the  J.  Fluid  Mech. ,  has 
been  completed. 


2.2  Stability  analysis 

A  paper  investigating  the  linear  stability  characteristics  of  the 
initial  region  of  a  two-dimensional  shear  layer,  including  the  effects 
of  the  splitter  plate  wake  as  well  as  the  effects  of  unequal  free  stream 
density,  was  presented  at  the  1987  AIAA  Aerospace  meeting  (Koochesfahani 
&  Frieler  1987).  The  important  point  of  this  paper  was  that  when  the 
effects  of  a  possible  velocity  profile  defect  (arising  from  the  splitter 
plate  wake)  are  taken  into  account,  it  is  found  that  the  flow  possesses 
two  unstable  modes.  In  particular,  in  addition  to  a  "shear”  mode  that 
can  be  recognized  as  evolving  from  the  classical  Kelvin-Helmholz 
instability,  a  second  "wake"  mode  is  also  present.  For  equal 
free-stream  densities,  the  wake  mode  is  not  as  unstable  as  the  shear 
mode.  For  large  values  of  the  density  ratio  and  for  the  case  where  the 
heavy  fluid  comprises  the  low  speed  stream  (p2/p^  >1),  however,  the 
wake  mode  can  be  considerably  more  unstable.  There  is  evidence  to 
suggest  that,  in  that  case,  the  importance  of  the  wake  mode  is  not 
confined  to  the  initial  region  but  extends  over  a  large  range  of  the 
flow.  This  paper  (Koochesfahani  &  Frieler  1987)  is  reproduced  as 
Appendix  A  of  this  report. 
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In  an  extension  of  this  work,  we  are  incorporating  the  effects  of 
compressibility  in  this  calculation.  Earlier  calculations 
(Gropengiesser  1970),  as  well  as  other  investigations  (Papamoschou 
1986),  suggest  that  such  effects  can  be  significant.  The  extension  of 
this  work  to  supersonic  stability  is  undertaken  as  part  of  the  on-going 
Ph.D.  research  of  Ms.  Mei  Zhuang. 


2.3  Penalty  ratio  effects 

Considerable  progress  has  been  made  in  our  experimental 
investigations  of  the  effects  of  free  stream  density  ratio  in  gas  phase, 
subsonic  shear  layers.  Specifically,  using  the  results  of  chemical 
reactions  at  both  high  and  low  stoichiometric  mixture  ratios,  i.e. 
4>  -  1  /8  and  *  8,  and  paying  particular  attention  to  avoid  finite 
Damkohier  number  effects,  we  have  been  able  to  determine  the  average 
high  speed  fluid  mixture  fraction  in  the  molecularly  mixed  fluid  (£m)  as 
a  function  of  free  stream  density  ratio,  for  density  ratios  p2/pi  in  the 
range  of  1/8  to  4  .  The  new  results  clarified  some  unresolved  questions 
that  were  raised  in  earlier  investigations  and  are  well  correlated  (in 
absolute  value)  with  estimates  computed  using  the  3imple  dependence  of 
the  entrainment  ratio  on  density  ratio  (Dimotakis  1986).  See  figure  1. 
We  are  investigating  whether  the  small  but  increasing  deviation  as  the 
density  ratio  increases  is  attributable  to  changes  in  the  large  scale 
flow  structure,  which  would  alter  the  entrainment  ratio,  or  increasing 
contributions  to  p(0,  the  probability  density  function  (PDF)  of  the 
mixed  fluid  fraction  £,  from  values  of  £  close  to  unity,  which  may  not 
have  been  adequately  represented  in  the  <p  -  1  /8  ,  8  "flip"  experiments. 
These  new  data  also  allow  us  to  estimate  the  dependence  of  the 
molecularly  mixed  fluid  mean  number  density  profile  across  the  layer,  as 
a  function  of  the  free  stream  density  ratio.  The  results  are  plotted  in 
figure  2.  The  important  conclusion  is  that  this  quantity  remains 
essentially  unaltered  as  the  free  stream  density  ratio  is  varied  over  a 
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wide  range.  Note  that  the  peak  values  in  this  calculation  are  slightly 
depressed  in  the  middle  of  the  layer  owing  to  the  reduction  of  the  total 
fluid  density  due  to  the  heat  release  and  need  to  be  properly  normalized 
with  the  corresponding  mean  number  density  profile  to  represent  the  mean 
mixed  fluid  mole  fraction  profile.  The  results  of  this  calculation  are 
depicted  in  figure  3. 

This  work  is  part  of  the  on-going  Ph.D.  research  effort  of  Mr. 
Cliff  Frieler. 


2.4  Dagkohler  number  effects 

The  gas  phase,  subsonic  low  heat  release  Damkohler  number 
experimental  investigation,  documented  in  the  GALCIT  Report  FM85-01  by 
Mungal  A  Frieler  (1985,  reproduced  also  as  appendix  B  in  Dimotakis, 

Broadwell  &  Leonard  1986),  has  been  accepted  for  publication  in  Comb.  & 

Flame. 

Pilot  calculations  of  Damkohler  number  effects  in  supersonic  shear  i 

layers  are  being  performed,  initially  undertaken  in  support  of  the 
design  effort  for  the  supersonic  shear  layer  facility.  We  subsequently 
decided  that  they  are  of  general  interest  in  their  own  right,  and  will 
be  presented  at  the  Joint  Propulsion  Meeting  in  San  Diego  this  June 
(Hall  A  Dimotakis  1987).  Briefly,  the  model  calculations  assume  as  a 
zeroth  order  flow/thermodynamic  model  that:  j 

1  .  the  chemical  kinetic  process  evolves  in  a  constant  pressure 

"balloon"  reactor,  whose  mass  is  increasing  linearly  with  time  | 

(a  consequence  of  the  entrained  reactant  flux  into  the  mixing  : 

layer),  ; 

2.  the  entrained  fluid  from  each  of  the  free  streams  is  mixed  and  * 

homogenizes  at  a  (single-valued)  composition  corresponding  to  9 

fj 
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i 

> 

the  entrainment  ratio, 

i  and, 

|  3-  that  it  is  entrained  in  a  thermodynamic  state  contributing  an 

|  energy  flux  from  each  of  the  free  streams  estimated  by 

!  computing  their  enthalpy  assuming  that  the  entrained  fluids  are 

|  brought  to  stagnation  conditions  in  the  convected  large  scale 

j  vortex  structure  frame. 

An  important  parameter  that  emerges  from  these  calculations  is  the 
I  dimensionless  ratio 


nAe 
nA(0 


which  measures  the  relative  importance  of  the  entrainment  flux  in  the 
chemical  kinetic  process.  In  this  expression,  nAe  is  the  lean  reactant 
species  entrainment  flux  (moles/sec),  tx  is  the  chemical  time  required 
to  consume  the  lean  reactant  in  a  constant  mass  reactor,  i.e.  in  the 
limit  of  small  and  nA(0)  is  the  initial  number  of  moles  of  the  lean 
reactant  species  in  the  reactor,  prior  to  the  initiation  of  the 
entrainment  process  at  t  »  0.  The  latter  quantity  is  intended  to 
represent  the  effective  amount  of  lean  reactant  that  is  rolled  up  in  the 
first  vortical  structure,  prior  to  the  origin  in  time  when  the  constant 
entrainment  flux  can  be  considered  to  have  begun. 


For  Qe  «  1  ,  the  chemical  kinetics  evolve  essentially  as  in  a 
constant  mass  reactor,  in  which  the  reaction  is  initiated  at  t  •  0  (with 
premixed  reactants)  and  proceeds  to  equilibrium  at  constant  pressure. 
For  »  1 ,  which  is  generally  the  fluid  mechanically  significant 
regime  for  high  Reynolds  number  shear  layers,  the  chemical  kinetics 
evolution  is  different,  reaching  an  asymptotic  behavior  that  is 
independent  of  Qg  for  large  values  of  that  parameter.  Sample 
calculations  for  the  H2/F2/N0  chemical  system  for  the  two  limiting 
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oases,  i.e.  «  i  (constant  mass  reactor;  solid  line),  and  »  1 
(entrainment-dominated  reactor;  dot-dashed  line),  are  depicted  in 
figure  4.  These  calculations  correspond  to  typical  subsonic  low  heat 
release  run  conditions  in  the  experiments  of  Mungal  &  Dimotakis  (1984) 
and  were  performed  using  the  CHEMKIN  software  library  package  (Kee  et  al 
’980)  and  the  Hp/f^/NC  chemical  reaction  system,  as  described  in  Mungal 
4  Frieler  '’985).  We  note  the  much  slower  approach  to  equilibrium  in 
the  entrainment-dominated  reactor;  a  consequence  of  the  continuous 
illation  of  tne  hot  product  mixture  by  the  relatively  cold  entrained 
fluids. 

Similar  calculations  were  performed  for  a  supersonic  shear  layer 

and  a  H2/F2/NO  reactant  system,  which  indicate  that  (qualitatively)  this 
system  behaves  quite  predictably.  In  particular,  as  the  entrained 
reactant  concentrations  and/or  enthalpy  are  raised,  the  effective 
kinetic  rate  increases  and  the  time  required  for  the  chemical  reaction 
to  be  driven  to  completion  decreases. 

The  H^/air  chemical  system  is  also  being  studied  in  similar  flow 
conditions  for  comparison  purposes,  as  well  as  for  its  significance  as  a 
potential  chemical  system  for  hypersonic  propulsion.  We  note,  however, 
that  for  this  system,  while  raising  the  static  temperature  (enthalpy)  of 
the  entrained  reactants  (preheating)  decreases  the  time  required  for 
ignition  (the  beginning  of  any  substantial  heat  release),  as  one  would 
expect,  it  can  also  increase  (possibly  substantially)  the  time  required 
for  completion  of  the  reaction.  This  is  a  consequence  of  the  fact  that, 
at  the  higher  absolute  temperatures  attained  under  these  conditions,  a 
large  fraction  of  the  enthalpy  in  that  system  can  be  tied  up  in 
relatively  small  concentrations  of  minor  but  energetic  species  that, 
especially  under  the  impeding  influence  of  the  cold  entrained  streams, 
reach  their  equilibrium  concentrations  rather  slowly. 

These  calculations  will  be  documented  in  the  Joint  Propulsion 
Meeting  AIAA  Paper  (Hall  ^  Dimotakis  1987). 
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2.5  Design  effort  for  the  supersonic  combustion  shear  layer  facility 

As  a  consequence  of  the  extension  of  the  scope  of  this  work  to 
variable  total  enthalpy  flows,  made  possible  by  co-funding  under  the 
AFOSR  F49629-86-C-01 1 3  URI  Contract,  the  overall  facility  design  had  to 
be  reviewed  with  this  important  new  specification  in  mind.  In 
particular,  the  design  specifications  for  several  critical  items  were 
modified.  This,  in  turn,  produced  revised  specifications  for  the 
high-pressure  gas  supply  vessel  and  the  fast-acting  valve/pressure 
regulator  (see  Dimotakis,  Broadwell  &  Leonard  1986,  section  2.5). 
Additionally,  two  design  concepts  for  the  test  section  are  being 
evaluated.  These  developments  are  briefly  described  below. 

A  final  specifications  document  has  been  written  for  the 
high-pressure  supply  vessel  (included  as  appendix  B).  This  document  has 
been  sent  to  local  manufacturers,  who  are  now  in  the  process  of 
formulating  official  contract  bids.  Following  receipt  of  the  bids  and 
on-site  visits  to  potential  vendors,  the  contract  will  be  awarded 
shortly  and  that  the  finished  tank  should  be  delivered  to  Caltech  by  the 
end  of  this  summer  or  early  fall. 

The  stringent  requirements  placed  on  the  pressure  regulator  (high 
pressure,  high  temperature,  large  flow  rate,  and  fast  response)  have 
eliminated  from  our  list  of  potential  vendors  all  the  commercial 
suppliers  of  conventional  pressure  regulators  known  to  us.  There  is, 
however,  reasonable  hope  that  the  Digital  Valve  (Digital  Valve  Corp.) 
will  essentially  meet  our  requirements .  This  valve  consists  of  an  array 
of  typically  12  valve  elements  of  different  sizes.  The  sizes  are 
arranged  in  a  binary  sequence.  Each  element  is  individually  operated 
from  a  digital  signal,  thus  in  principle,  within  the  resolution  of  one 
part  per  212  (4,096),  any  desired  total-area-versus-time  curve  can  be 
programmed  into  the  valve.  The  valve  model  currently  under 
consideration  (2"-DIGIC0N-RC)  meets  the  pressure  specifications.  The 
temperature  specifications  are  being  discussed  with  the  manufacturer  to 
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explore  variations  in  the  basic  design  permitting  the  accommodation  of 
our  temperature  specifications1-.  We  estimate  that  two  such  valves, 
operated  in  parallel  (also  increasing  the  resolution  thereby),  will  have 
adequate  flow  capacity.  There  are  additional  questions,  however, 
concerning  the  temporal  resolution  capabilities  of  the  valve.  The 
quoted  response  time  of  75  msec  for  the  opening/closing  of  each  element 
made  us  apprehensive  about  the  smoothness  of  the  total-area-versus-time 
output  and  the  resulting  flow  quality.  Consequently,  an  analytical 
model  of  the  valve  dynamics  was  devised,  in  which  the  area  change  of  a 
single  element  during  the  opening/closing  cycle  was  represented  by  a 
smooth  curve.  Our  analysis  predicts  that  a  very  smooth  output  can  be 
obtained  when  the  valve  configuration  is  updated  at  time  intervals  close 
to  the  response  time.  See  figure  5.  In  summary,  the  prospects  of  using 
the  Digital  Valve  are  promising,  provided  the  manufacturer  can  overcome 
the  temperature  limitations  in  the  present  design  and  that  he  (and  we) 
confirm  our  timing  predictions. 

Finally,  we  are  considering  two  candidate  nozzle  configurations  for 
producing  a  variable  exit  Mach  number  in  the  range  from  1.5  to  3.5  .  In 
particular,  we  are  evaluating  the  relative  merits  of  a  symmetric  2-D 
nozzle  with  replaceable  contours,  or  an  asymmetric  2-D  nozzle  with  a 
sliding  block.  The  former  has  the  advantage  of  design  simplicity  and 
the  disadvantage  of  the  need  to  fabricate  a  nozzle  with  different 
contours  for  each  different  Mach  number  and  gas  (Y)  combination.  The 
latter  offers  the  flexibility  of  a  continuous  Mach  number  range,  with 
some  (minor)  compromises  in  flow  quality,  however,  as  well  as  some 
complications  in  the  over-all  test-section  design  and  producibility. 
Additionally,  the  asymmetric  nozzle  needs  to  be  about  twice  as  long  as 
the  symmetric  one,  resulting  in  a  splitter  plate  trailing-edge  (high 
speed)  boundary  layer  momentum  thickness  (9^  that  will  be  roughly  twice 
as  thick.  This  is  an  important  consideration,  as  the  flow  streamwise 


As  of  this  writing,  the  high  temperature  (600°F)  specification  is 
causing  some  difficulties,  as  the  stock  item  is  fabricated  with  soft 
seals  in  the  air  cylinder  actuators. 
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run  (L)  is  more  or  less  fixed  by  other  constraints  and  we  would  like  to 
operate  at  values  of  L/91  that  are  as  large  as  possible  and,  in 
particular,  no  less  than  10^,  or  so.  We  are  presently  acquiring  and 
adapting  computer  flow  codes  to  assist  us  in  the  design  and  relative 
evaluation  of  these  nozzles.  A  decision  as  to  which  configuration  will 
be  implemented  is  expected  in  the  next  few  months. 

This  work  is  part  of  the  on-going  Ph.D.  research  effort  of  Mr. 
Jeffery  Hall,  in  collaboration  with  Dr.  Dimitri  Papamoschou  and  Mr. 
Cliff  Frieler. 


2.6  Liquid  phaae  reacting  and  non-reacting  shear  layers 


The  archival  documentation  of  our  earlier  efforts  in  this  area  has 
just  appeared  in  the  J.  Fluid  Mech.  (Koochesfahani  &  Dimotakis  1986). 


.  . ...  -  -  - 


2.7  Mew  shear  layer  al»lng  model 

The  approach  adopted  in  the  new  model  is  that  of  viewing  an 

Eulerian  slice  of  the  spatially  growing  shear  layer,  at  a  downstream 

station  at  x,  and  imagining  the  instantaneous  interface  between  the  two 
interdiffusing  and  chemically  reacting  fluids  as  well  as  the  associated 
strain  field  Imposed  on  that  interface.  It  is  recognized  that  both  the 
Eulerian  state  and  the  local  behavior  of  that  interface  are  the 

consequence  of  the  Lagrangian  shear  layer  dynamics  from  all  relevant 
points  upstream  of  the  station  of  interest  at  x.  It  is  assumed, 

however,  that  this  upstream  history  acts  in  such  a  manner  as  to  produce 
a  self-similar  state  at  x,  whose  statistics  can  be  described  in  terms  of 
the  local  parameters  of  the  flow.  In  particular,  it  is  assumed  that  a 
Kolmogorov  cascade  process  has  been  the  appropriate  description  of  the 
upstream  dynamics,  leading  to  the  local  Eulerian  spectrum  of  scales  and 
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associated  strain  rate  field  at  x. 

The  justification  for  this  approach  is  that  while  the  large  scale 
dynamics  are  all  important  in  determining  such  things  as  the  growth  rate 
and  entrainment  ratio  into  the  spatially  growing  shear  layer,  the 
predominant  fraction  of  the  interfacial  area  is  associated  with  the 
smallest  scales,  which  can  perhaps  be  adequately  dealt  with  in  terms  of 
universal  similarity  laws.  The  large  scales,  therefore,  are  viewed  as 
feeding  the  reactants  that  are  entrained  at  some  upstream  station  into 
the  smaller  scale  turbulence  at  the  appropriate  rate.  These  reactants 
subsequently  get  processed  by  the  evolution  of  the  cascade  processes 
upstream  to  produce  the  local  spectrum  of  scales  at  x  (see  discussion  in 
Broadwell  &  Dimotakis  1986).  Thi3  conceptual  basis  is  also  aided  by  the 
notion  of  a  conserved  scalar,  according  to  which  the  state  of  diffusion 
and  the  progress  of  an  associated  chemical  reaction,  in  the  limit  of 
fast  (diffusion-limited)  chemical  kinetics,  is  completely  determinable 
by  the  local  ( Euler ian)  state  of  the  conserved  scalar  (3ee,  for  example, 

Bilger  1980). 

An  important  part  of  the  proposed  procedure  is  the  normalization 
that  is  imposed  on  the  statistical  weight  (contribution)  of  each  scale  X 
to  the  total  amount  of  molecularly  mixed  fluid  and  associated  chemical 
product.  This  is  done  via  the  expected  interfacial  surface  per  unit 
volume  ratio  that  must  be  assigned  to  each  scale  X.  When  totalled  over 
all  scales,  these  statistical  weights  must  add  up  to  unity. 

The  results  are  first  obtained  conditional  on  a  uniform  value  of 
the  dissipation  rate  e  .  An  attempt  to  incorporate  and  assess  the 
effects  of  the  fluctuations  in  the  local  dissipation  rate,  i.e. 
e(x,t)  ,  is  made  by  folding  the  conditional  results  over  a  probability 
density  function  for  e  . 


In  a  similar  vein,  a  refinement  of  the  entrainment  ratio  idea  is 
incorporated,  in  which  it  is  recognized  that  the  large  scale  spacing  l/x 
is  a  random  variable  and  that  therefore,  by  the  force  of  the  arguments 
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in  Dimotakis  (1986),  the  entrainment  ratio  must  also  be  treated  as  a 
random  variable  of  the  flow.  Accordingly,  the  results  are  obtained 
conditional  on  a  given  value  of  the  entrainment  ratio  E,  and  are 
subsequently  folded  over  the  expected  distribution  of  values  of  E  about 
its  average  value  E. 

In  the  calculations,  it  is  assumed  that  the  molecular  diffusi vities 
for  all  relevant  species  are  equal  to  each  other,  but  not  necessarily 
equal  to  the  kinematic  viscosity.  Heat  release  effects  and  temperature 
dependence  effects  of  the  molecular  transport  coefficients  are  also 
ignored.  This  is  appropriate  for  the  liquid  phase  measurements  of 
Koochesfahani  A  Dimotakis  (1986),  and  may  be  adequate  for  the 
description  of  the  gas  phase  measurements  of  Mungal  A  Dimotakis  (1984) 
and  the  Reynolds  number  study  of  Mungal  et  al  (1985).  The  issue  of  heat 
release  effects  on  the  flow  was  specifically  addressed  elsewhere  (see 
Hermanson  et  al  1987).  In  computing  the  temperature  corresponding  to 
the  heat  released  in  the  reaction,  equal  heat  capacities  are  also 
assumed  for  the  two  fluids  brought  together  within  the  mixing  zone. 
While  some  of  these  assumptions  were  not  necessary,  they  allowed 
calculations  to  be  performed  in  closed  form  permitting,  in  turn,  the 
examination  of  the  dependence  of  the  results  on  the  various 
dimensionless  parameters  of  the  problem.  The  proposed  procedure  assumes 
that  the  relevant  statistics  of  the  velocity  field  are  known  (or  can  be 
estimated)  and  computes  the  behavior  of  the  passive  scalar  process  in 
response  to  that  velocity  field.  Finally,  the  procedure  is  "closed"  in 
that  it  yields  the  (absolute)  chemical  product  volume  fraction  6p/s  in 
the  shear  layer  at  x,  with  no  adjustable  parameters. 

The  prediction  for  the  chemical  product  thickness  6pj ($)/&,  where 
<5pi  is  given  by 


( $) 
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where  Cp(y,$)  is  the  chemical  product  concentration,  cqi  is  the  reactant 
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concentration  in  the  high  speed  stream,  6  is  the  total  (1$)  thickness  of 
the  layer  and  <p  is  the  stoichiometric  mixture  ratio  are  depicted  in 
figure  6  .  The  data  (circles)  are  from  the  gas  phase  measurements  of 
Mungal  A  Dimotakis  (1984).  The  correct  prediction  (in  absolute  value) 
is  noteworthy. 

The  predictions  for  the  chemical  product  mean  volume  fraction, 
<5p/s ,  where 


6p 


f  AT(y,f) 
L  ATflB(f) 


(3) 


AT(y,$)  is  the  mean  temperature  rise  profile,  and  where  ATflm($)  is  the 
adiabatic  flame  temperature  at  the  stoichiometric  mixture  ratio  <p,  are 
plotted  (solid  line)  in  figure  7  versus  the  stoichiometric  mixture 
fraction,  i.e. 


«♦ 
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for  the  Mungal  &  Dimotakis  (1984)  data  depicted  in  figure  6  (note  that 
SP1  ”  £$6p).  The  corresponding  predictions  are  also  plotted  (dashed 
line)  for  the  liquid  phase  point  of  Koochesfahani  &  Dimotakis  (1986),  at 
♦  -  10  and  a  comparable  Reynolds  number.  Note  that  the  Schmidt  number 
dependence  of  the  chemical  product  volume  fraction  6p/6  is  also 
predicted  correctly. 


Finally,  the  prediction  for  the  Reynolds  number  dependence  of  5p/j 
is  depicted  in  figure  8  and  compared  with  the  gas  phase  measurements  of 
Mungal  et  al  (1985)  and  the  liquid  phase  measurements  of  Koochesfahani  A 
Dimotakis  (1986). 

Some  of  the  features  of  this  model  and  its  predictions,  which 
distinguish  it  from  the  model  proposed  by  Broadwell  &  Breidenthal 
(1982),  and  as  updated  in  the  discussion  in  Broadwell  A  Mungal  (1986), 
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are  listed  below,  are: 

1.  There  are  no  adjustable  parameters  in  this  model. 

2.  5p/5,  the  volume  fraction  in  the  shear  layer  region  occupied  by 
chemical  product  tends  to  zero  as  the  Schmidt  number  tends  to 
infinity  at  fixed  Reynolds  number  (there  is  no  explicit 
procedure  for  computing  this  dependence  in  the 
Broadwell-Breidenthal  model  for  large  Schmidt  numbers). 

3.  <Sp/6  is  also  predicted  to  tend  to  zero  in  the  limit  of  infinite 

Reynolds  numbers,  albeit  slowly  (logarithmically).  The 

corresponding  limit  in  the  Broadwell-Breidenthal  model  is 
finite  (see  discussion  in  section  3-5). 

I 

I 

>  This  work  will  be  presented  at  the  U.S.  -  France  Joint  Workshop  on 

I 

[  Turbulent  Reacting  Flows  this  July  (Rouen,  France)  and  will  be  issued  in 

the  interim  as  a  GALCIT  report  (Dimotakis  1987). 

3.0  TURBULENT  JET  MIXING  &  CHEMICAL  REACTIONS 

Parts  of  this  work  have  been  cofunded  by  the  Gas  Research 
Institute,  the  Office  of  Naval  Research  and  the  National  Bureau  of 
Standards. 


3.1  Gaa  phase  Mixing  at  low  to  Moderate  Reynolds  nuMbers 

To  increase  the  understanding  of  the  mixing  processes  in  gas  phase 
turbulent  jets,  an  investigation  of  the  entire  range  of  concentration 
length  and  time  scales  has  been  undertaken  at  low  to  moderate  Reynolds 
numbers.  In  the  first  year  of  this  effort  a  laser  Rayleigh  scattering 
diagnostic  was  developed  and  the  design  of  an  experiment  to  exploit  the 
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power  of  this  technique  was  begun.  In  the  second  year,  the  design  was 
completed,  the  experiment  was  constructed  and  performance  testing  was 
begun.  Last  year  was  devoted  to  the  resolution  of  problems  identified 
by  the  performance  tests  and  the  improvement  of  some  aspects  of  the 
experiment  with  new  technologies.  A  few  new  runs  were  made  at  the  full 
capabilities  of  the  apparatus  with  some  interesting  results.  A  method 
of  flow  visualization  and  some  new  data  processing  software  have  also 
been  developed. 

In  last  year's  report  (Dimotakis,  Broadwell  &  Leonard  1986,  section 
4.1),  it  was  stated  that  the  coflow  became  unsteady  after  about  15 
seconds.  This  problem  was  attributed  to  a  cooling  trend  in  the  coflow 
gas  in  the  course  of  a  run.  Because  the  coflow  gas  enters  at  the  top  of 
the  experiment,  a  cooler  (denser)  layer  of  gas  was  able  to  form  over  a 
warmer  (lighter)  gas.  This  configuration  is  gravitationally  unstable  to 
the  formation  of  buoyancy  currents.  A  feedback  temperature  control 
system  was  designed,  fabricated  and  installed  to  address  this  problem. 
Nevertheless,  the  unsteady  buoyancy  currents  in  the  test  section  were 
not  entirely  eliminated.  The  remaining  problem  was  diagnosed  to  be  a 
consequence  of  small  isothermal  density  gradients,  which  could  be 
established  in  the  facility  prior  to  the  initiation  of  the  run.  A  new 
purging  procedure  was  subsequently  adopted  which  solved  this  last 
problem.  The  coflow  is  now  stable  for  approximately  one  minute,  which 
is  long  enough  to  accommodate  the  scheduled  runs. 

Preliminary  calculations  had  indicated  that,  with  the  advent  of  a 
new  line  of  high  quality  photodiodes,  the  possibility  of  replacing  the 
photomultiplier  tube  (PMT)  as  the  photodetector  should  be  considered. 
The  photodiode  has  the  advantage  of  higher  quantum  efficiency  (0.6  to 
0.7  in  the  wavelength  range  of  interest),  relative  to  that  available 
from  the  PMT  (typically  0.1).  The  disadvantage  of  a  photodiode  is  that 
it  does  not  amplify  its  signal  (Gain  -  1),  as  does  the  PMT  (Gain  -  10^), 
and  therefore  must  be  mated  to  a  very  high  quality,  low  noise 
( trans impedance )  amplifier  if  it  is  to  compete.  A  typical  signal 
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current  for  this  experiment,  for  example,  would  be  about  0.3*10 A  from 
a  photodiode  and  about  0.3»10_^A  from  the  rMT.  The  design  of  the 
matched  trans impedance  amplifier  was  completed  by  Dr.  Daniel  Lang  and  a 
performance  calculation  was  undertaken  for  the  proposed  matched 
photodiode/amplifier  system.  The  results  were  promising  enough  to 
suggest  that  building  and  testing  a  prototype  was  warranted.  Following 
a  few  modifications  to  decrease  stray  capacitance  and  increase 
bandwidth,  the  photodiode  system  performance  was  indeed  proven  better 
than  the  corresponding  PMT  system  performance.  The  photodiode  system  is 
also  a  more  attractive  detector  because  it  is  inexpensive,  mechanically 
tough  and  optically  rugged.  It  is  driven  by  a  low  voltage  (+15  Volts), 
compared  to  the  PMT  (-  800  Volts),  which  results  in  a  safety  advantage 
since  possible  sparks  are  of  some  concern  when  combustible  gases  are 
present  in  the  measurement  environment  (most  of  the  Jet  gases  used  in 
this  study  are  hydrocarbons).  These  facts  made  the  decision  to  replace 
the  PMT  system  with  that  using  the  photodiode  compelling.  Figure  9  is  a 
comparison  plot  of  the  measured  power  spectrum  from  each  of  the  systems 
when  illuminated  by  the  Rayleigh  scattered  light  from  argon  gas  at  1 
atmosphere  pressure.  The  scattered  light  was  collected  from  a  200  pm 
length  of  a  collimated  15  Watt  Ar+  laser  beam  in  both  cases.  As  can  be 
seen  from  these  data,  the  spectrum  below  about  30  Hz,  for  both  systems, 
is  dominated  by  1/f  noise  (from  the  laser).  The  main  result  to  note  is 
that  from  about  30  Hz  to  about  20  kHz,  the  bandwidth  of  interest  for 
these  studies,  the  photodiode  system  noise  spectrum  falls  below  that  of 
the  PMT  system. 

An  effort  was  also  made  to  update  the  existing  data  acquisition 
software  to  take  full  advantage  of  the  increased  A/D  conversion 
capability,  computing  power,  and  memory  available  on  the  laboratory 
computer  used  for  this  experiment.  In  particular  this  allowed  the 
length  of  a  data  set  for  each  run  to  be  comfortably  increased  by  a 
factor  of  2,  to  2^  (-  500,000)  individual  concentration  measurements. 
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For  this  investigation,  the  power  spectrum  of  concentration 
fluctuations  is  an  important  diagnostic  of  the  turbulent  state  of  jet 
mixing.  In  order  to  calculate  a  spectrum  over  the  possible  range  of 
frequencies  it  was  necessary  to  develop  algorithms  to  calculate  the  high 
and  low  frequency  parts  of  the  spectrum  separately.  These  routines  were 
tested  with  some  preliminary  data  and  then  used  to  compute  the  noise 
spectra  (figure  9)  mentioned  above. 

To  ensure  that  the  experimental  apparatus  is  producing  a  "standard" 
jet,  it  is  necessary  to  check  that  the  half-angle  of  the  cone  which 
contains  the  jet  fluid  falls  within  the  range  of  accepted  values 
01°  i  a  t  12.5°),  i.e. 


6(x) 
2  x 


tan  a 


where  S(x)  is  the  local  jet  diameter  at  the  edge  of  the  jet  fluid 
containing  region.  This  angle  can  be  estimated  experimentally  from 
photographs  of  the  jet  flow.  Because  the  gases  used  in  this  study  are 
clear,  shadowgraph  photography  was  the  obvious  choice.  In  this  case, 
the  shadowgraphy  relies  on  the  same  difference  in  the  refractive  index 


of  the  jet  and  reservoir  gases  exploited  in  the  Rayleigh  scattering 
measurements.  A  sample  photograph  is  included  as  figure  10.  The 
spreading  half-angle  for  this  jet  is  about  11.5°.  The  field  of  view  is 
from  x/d  -  0  (top  of  photograph)  to  x/d  -  40  (bottom  of  photograph)  jet 
exit  diameters.  The  jet  gas  is  propylene,  the  reservoir  gas  is  argon 
and  the  Reynolds  number  is  5,000.  The  slight  non-uniformity  in  the 
background  light  level  arises  because  the  arc-lamp  light  source  used 
does  not  produce  a  perfect  spherically  diverging  beam. 

With  the  coflow  problems  solved  and  the  entire  data  acquisition 
system  redesigned,  reprogramed  or  upgraded,  the  high  quality 
measurements  that  are  focus  of  this  work  were  begun.  A  power  spectrum 
of  the  concentration  fluctuations  in  the  jet  at  x/d  »  40  has  been 
Included  as  figure  11.  The  diagonal  line  has  a  slope  of  ~5/3*  The 
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frequency  axis  is  normalized  by  the  estimated  local  large  scale  time, 
estimated  here  by  means  of  the  expression 


5 

U 


where  U  -  U(x)  is  the  local  mean  jet  velocity  at  the  centerline.  This 
is  estimated  to  be  about  0.6  seconds  at  this  measuring  station.  A 
comparison  with  the  noise  spectra  from  the  sensors  shows  that  the 
concentration  fluctuations  in  the  jet  are  3  to  4  orders  of  magnitude 
larger  that  the  system  noise.  Recognizing  that  the  turbulent  cascade 
may  be  underdeveloped  for  this  Reynolds  number,  a  second  and  potentially 
more  important  observation  is  that  the  estimated  Kolmogorov  scale  mean 
passage  frequency,  i.e. 


'k  -  f 

aK 


where  the  Kolmogorov  scale  XK  is  here  estimated  by, 
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and  the  observed  "end"  of  the  -5/3  spectral  slope  differ  by  almost  a 
factor  of  23  (at  the  Reynolds  number  of  5,000  for  this  run).  This  means 
that  the  smallest  scalar  diffusion  scales  (Sc  »  1  here)  are  much  larger 
than  our  design  estimates.  We  consequently  decided,  so  as  to  exploit 
the  full  capabilities  of  our  diagnostics,  to  reconfigure  the  experiment 
so  as  to  permit  higher  Reynolds  numbers  to  be  attained.  This  will 
involve  increasing  the  volume  flux  from  the  coflow  and  designing  and 
fabricating  a  new  jet  nozzle  with  a  smaller  exit  diameter.  The  new 
nozzle  should  allow  the  jet  Reynolds  number  to  be  increased  to  almost 
25,000  while  observing  the  strict  spatial  and  temporal  resolution 
constraints  of  the  experiment. 
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We  plan  to  complete  the  measurements  and  data  analysis  at  the 
current  Reynolds  number  (5,000),  and  then  move  to  the  new  configuration 
to  explore  Reynolds  number  effects  in  the  coming  year. 

This  work  is  part  of  the  on-going  Ph.D.  research  effort  of  Mr. 
David  Dowling. 


3-2  High  pressure  jet  combustion  facility 

The  High  Pressure  Reactant  Vessel  is  designed  for  studies  of  high 
Reynolds  number  turbulent  reacting  jets,  at  both  reduced  and  elevated 
pressures.  In  addition  to  the  capability  of  using  conventional  gaseous 
fuels,  the  vessel  has  been  designed  to  handle  more  exotic  reactants  like 
fluorine,  nitric  oxide,  and  hydrogen.  This  capability  will  allow  a 
systematic  study  of  the  effects  of 

1.  stoichiometric  mixture  ratio, 

2.  Jet/reservoir  fluid  density  ratio, 

3.  heat  release, 

4.  finite  kinetic  rate  (Damkohler  number), 
as  well  as 

5.  buoyancy, 

over  a  wide  range  of  Reynolds  numbers. 

The  High  Pressure  Reactant  Vessel  (HPRV)  fabrication  was  completed 
and  the  vessel  was  assembled  in  the  lab  where  it  was  helium  leak  tested 
and  vacuum  tested  to  our  satisfaction.  Figure  12  shows  the  vessel 
installed  in  the  lab. 
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The  gas  delivery  system  will  be  modeled  after  the  H2/f2  combustion 
shear  layer  delivery  system.  In  order  to  ensure  the  safety  of  this 
design  in  the  new  facility,  which  will  see  much  higher  concentrations  of 
fluorine  at  higher  pressures  than  the  shear  layer  facility,  a  materials 
test  was  performed  at  the  TRW  Chemical  Laser  Lab  in  Redondo  Beach.  The 
materials  to  be  exposed  to  fluorine  in  the  new  lab  were  subjected  to  5 
times  the  highest  concentration  of  fluorine  we  expect  to  see  in  the  lab 
at  the  highest  pressure.  All  materials  survived  the  tests,  indicating 
that  we  should  proceed  with  our  gas  delivery  design  as  planned.  Most  of 
the  components  of  this  delivery  system  have  been  purchased  and 
construction  of  the  control  panel,  schematically  shown  in  figure  13, 
will  begin  soon. 

Different  diagnostic  tools  are  being  considered  for  the  initial 
nonreacting  cold  runs  and  the  ensuing  reacting  hot  runs.  Shadowgraphy 
flow  visualization  will  be  a  first  diagnostic.  For  the  hot  runs,  the 
main  (additional)  diagnostic  for  the  initial  experiments  will  be  cold 
wire  thermometry,  similar  to  the  ones  currently  in  use  in  the  H2/f2 
shear  layer  facility.  Preliminary  tests  have  also  been  conducted  by  Mr. 
Cliff  Frieler  using  more  sensitive  wires,  which  are  also  character ized 
by  considerably  improved  time/space  resolution,  of  smaller  diameter 
(0.6  pm,  versus  the  2.5  pm  diameter  wires  presently  in  use  in  the  shear 
layer  experiments).  The  smaller  wires  are  being  fabricated  by  Mr.  Earl 
Dahl.  The  amplifier  circuits  for  these  probes  have  been  constructed  and 
will  allow  16  channels  of  high  temporal/spatial  resolution  temperature 
to  be  measured.  We  hope  to  begin  cold  runs  by  the  end  of  the  summer. 
Subsequent  experiments  will  also  benefit  from  the  diagnostics  and 
instrumentation  developments  presently  in  progress  under  the  guidance  of 
Dr.  R.  Miake-Lye  and  Dr.  D.  Lang  and  sponsored,  in  part,  by  the 
AFOSR-85~01 53  URI  contract. 

The  design  of  this  facility  begun  under  the  guidance  of  Dr.  G. 
Mungal  (presently  at  Stanford),  primarily  under  funding  from  the  Gas 
Research  Institute.  The  work  described  here  is  part  of  the  on-going 
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Ph.D.  research  effort  of  Mr.  R.  Gilbrech. 


3-3  Laser  scattering  exper Intents  In  hydrocarbon  flames 

A  paper  was  recently  published  in  Combustion  and  Flame  (Miake-Lye  3t 
Toner  1987,  included  as  appendix  C  in  this  report)  which  describes  our 
experimental  studies  of  large  buoyant  diffusion  flame  using  laser  soot 
scattering.  In  that  paper,  arguments  were  presented  justifying  the  use 
of  soot  as  a  marKer  for  the  combustion  interface  for  the  flame  under 
study.  In  particular,  we  argue  that  the  soot  is  localized  in  a  thin 
layer  close  to  and  locally  following  the  reaction  interface  surface. 

Based  on  our  r-t  digital  image  data,  we  have  calculated  the  soot 
intermittency  (-  soot  volume  fraction),  or  fraction  of  time  that  soot  is 
present,  at  each  radial  station.  By  determining  the  instantaneous  flame 
boundary,  we  have  also  calculated  tn..  conditional  soot  intermittency, 
which  we  define  to  be  the  fraction  of  time  that  soot  is  present  within 
the  flame  boundary  at  any  radial  station.  The  important  result  is  the 
observation  that  the  conditional  soot  intermittency  is  nearly  constant 
across  the  flame.  These  results  lead  us  to  a  picture  of  the  flame  made 
up  of  thin  laminar  diffusion  flamelets  distributed  fairly  uniformly 
within  the  puff-like  large  scale  envelope  of  the  flame.  The  flamelets 
are  continuous  laminar  diffusion  flames  near  the  burner,  but  become 
increasingly  convoluted  and  broken  (non-contiguous )  as  the  flame  tip  is 
approached. 

Subsequent  soot  scattering  experiments  have  concentrated  on  axially 
directed  imaging  in  the  same  large  buoyant  diffusion  flame.  These 
experiments,  which  were  proposed  in  the  last  annual  report,  were  carried 
out  at  various  fuel  flow  rates  with  the  beam  directed  parallel  to  the 
axi3  of  the  burner  at  each  of  three  radial  locations.  The  puffing 
periodicity  is  evident  in  the  reconstructed  z-t  images  and  the 
celerities  of  these  structures  can  be  measured  from  curve  fits  to  the 
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processed  digital  image  data. 

Data  analysis  is  in  progress.  Some  of  the  software  developed  for 
the  radial  data  analysis  has  been  modified  for  the  present 

investigation,  new  programs  have  been  written  (due  to  the  different 
structure  of  the  axial  images),  and  software  has  also  been  written  for 
use  on  a  new  image  processing  system.  Work  to  date  ha3  begun  to 
elucidate  the  accelerating  flow  region  close  to  the  burner,  driven  by 
the  heat  release  and  marked  by  the  concomitant  soot  sheets.  Not 
surprisingly,  the  structures  persist  longer  in  the  z-t  diagrams  for 

higher  fuel  flow  rates  (consistent  with  longer  flame  lengths)  and, 
further,  the  celerity  of  a  structure  at  a  given  axial  station  depends  on 
the  fuel  flow  rate.  Celerities  are  measured  to  range  from  50  cm/s  at 

z/D  -  0.2  to  190  cm/s  at  z/D  -  1  and  beyond.  Analysis  should  be 

completed  soon  and  a  report  of  these  results  will  be  written  sometime 
during  the  current  contract  year. 

In  support  of  our  ongoing  efforts  to  understand  flame  structure,  a 
feasibility  experiment  was  performed  using  an  additional  flow 
visualization  technique.  This  technique  makes  use  of  titanium 
tetrachloride  (TiCljj)  which,  in  gas  phase,  reacts  with  water  vapor  to 
form  Ti02  particles  (e.g.  Roqueraore  et  al.  1986),  which  scatter  laser 
light  strongly.  The  particles  will  be  present  anywhere  in  the  flow 

where  TiCl^  and  water  vapor  have  mixed,  so  by  adding  TiCl^  to  the  fuel 

or  ambient  air  and  water  vapor  to  the  same  fluid  or  the  other  fluid,  a 

variety  of  flow  labeling  schemes  can  be  achieved.  The  results  of  such 

experiments  can  complement  the  information  derived  from  soot  scattering 
experiments. 

The  feasibility  study  was  carried  out  with  the  50  cm  burner 
apparatus  by  including  a  partially  filled  TiCl^  container  in  the  fuel 
supply  line.  The  fuel  passed  over  the  surface  of  the  liquid  TiCl^, 
which  introduced  (an  as  yet  unquantified  amount  of)  TiC 1 ^  vapor  into  the 
fuel  gas.  When  the  TiCljj  vapor  encountered  water  vapor,  either  due  to 
unreacting  mixing  with  ambient  air  or  upon  combustion  of  the  fuel,  Ti02 
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particles  were  formed.  In  a  diffusion  flame,  soot  particles  are  also 
often  formed  at  the  reaction  interface,  so  the  feasibility  study  wa3 
intended  to  determine  whether  the  two  sources  of  particles  could  be 
distinguished. 

The  burner  was  lit  with  a  fuel  flow  rate  corresponding  to  -  1 90  kW 
heat  release,  a  value  (the  highest)  used  in  our  previous  studies.  The 
flame  was  run  with  and  without  the  TiCl^  seeding  and  an  increase  in  the 
luminosity  by  a  factor  of  3  -  5  was  apparent  with  TiCltj.  This  was 
determined  by  noting  the  attenuation  required  to  equalize  exposures  with 
and  without  the  TiCl^  seeding.  The  flame  luminosity  also  changed  color 
from  the  typical  yellow  color  due  to  hot  soot  to  an  intense  white  color. 

This  feasibility  study  leads  us  to  be  cautiously  optimistic  about 
the  prospects  for  using  TiCly  as  a  flow  visualization  tool  in  our 
continuing  flame  structure  experiments.  It  appears  that  when  it  is 
added  to  the  fuel  in  the  manner  described  above,  there  is  a  sufficiently 
large  increase  in  the  flame  luminosity  that  the  luminosity  due  to  the 
soot  will  be  a  fraction  (20  -  30%)  of  the  total.  This  should  permit  the 
type  of  selective  flow  labeling  that  we  envision  as  required  for  our 


ongoing  experiments. 


These  studies  were  undertaken  by  Dr. 

collaboration  with  Dr.  S.  Toner. 


Miake-Lye 


3. A  Liquid  phase  jet  mixing  at  small  scales 

Work  is  in  progress  to  study  the  mixing  interface  and  small-scale 
structure  of  turbulence  in  a  liquid  phase,  axisymmetric ,  turbulent  jet. 
The  diagnostic  method  used  is  Laser  Induced  Fluorescence,  which  has 
several  desirable  attributes  including  the  potential  (for  liquid  phase 
fluorescence)  of  very  high  signal-to-noise  ratio,  high  spatial 
resolution,  and  an  intrinsic  connection  with  the  mixing/diffusion 
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process.  Analysis  of  preliminary  data  has  been  undertaken  using  fractal 
and  spectral  methods. 

The  fractal  processing  algorithms  developed  are  based  on  the 
"box-counting'1  scheme,  which  we  have  implemented  in  both  one  and  two 
dimensions.  Issues  addressed  in  the  fractal  algorithms  have  included 
the  effect  of  the  thresholding  process  (which  determines  the  boundary  of 
the  region  to  be  "dimensioned"),  the  impact  of  finite  record  length, 
comparison  of  1 -D  vs  2-D  algorithms,  and  whether  the  turbulent  jet 
interface  is  in  fact  describable  by  the  notion  of  a  characteristic 
fractal  dimension.  Results  are  relatively  insensitive  to  the 
thresholding  level  and  on  whether  one  or  two  dimensional  algorithms  are 


It  was  found  that,  at  least  within  the  range  of  (the  small)  scales 
examined,  the  Jet  does  not  seem  to  possess  a  characteristic  fractal 
dimension.  Additionally,  end  effects  from  finite  record  lengths  can 
produce  substantial  noise  for  short  records  (i.e.  128  elements).  The 
smallest  scale  examined  is  estimated  to  be  a  fraction  of  the  Kolmogorov 
scale,  while  the  width  of  the  records  was  approximately  a  tenth  of  the 
jet  diameter  (consisting  of  512  elements).  While  the  scales  used  may  be 
either  too  small  (an  unusual  difficulty)  or  their  range  too  narrow,  and 
the  Reynolds  number  (1,000  -  3,000)  may  perhaps  be  rather  low,  the  lack 
of  a  characteristic  fractal  dimension  is  in  conflict  with  the  only  other 
such  account  in  the  literature  (Sreenivasan  &  Meneveau  1986),  where  the 
difficulties  of  end  effects  are  not  addressed  at  all. 

These  preliminary  investigations  have  prompted  us  to  consider 
additional  experiments,  tailored  to  the  collection  of  data  specifically 
suited  for  these  techniques.  The  new  experiments  will  be  optimized  for 
higher  Reynolds  numbers  (-  20,000),  which  should  be  sufficiently  large 
for  more  definitive  conclusions,  and  with  a  design  spatial/temporal 
resolution  of  the  order  of  a  few  Batchelor  (scalar  diffusion)  scales,  a 
region  never  Investigated  in  a  laboratory  jet  at  high  Schmidt  number,  to 
our  knowledge.  We  are  hoping  that  the  results  from  these  experiments 
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will  complement  in  many  important  ways  the  gas  phase  mixing  measurements 
under  way  (section  3.1)  and  help  clarify  important  issues  in  our 
understanding  of  turbulent  jet  mixing  phenomena. 


This  work  is  part  of  the  on-going  Ph.D.  research  of  Mr.  Paul 


Miller. 


3.5  Turbulent  jet  Modeling  efforts 

During  the  past  year,  modeling  efforts  based  on  the  ideas  described 
in  Broadwell  &  Breidenthal  (1982)  have  focused  in  two  areas: 

1 .  attempts  to  clarify  the  connection  of  the  model  foundations  to 
the  classical  ideas  of  Kolmogorov 


2.  the  development,  in  quantitative  and  useable  form,  of  an 
axi-symmetric  jet  model  corresponding  to  that  for  the  shear 
layer  (Broadwell  &  Mungal  1986). 

The  first  task  has  centered  around  the  connection  between  kinetic 
energy  dissipation  and  scalar  mixing.  While  no  definite  conclusions 
have  yet  been  reached,  some  progress  can  be  reported.  The  calculation 
by  Broadwell  A  Mungal  (1986),  of  the  fractional  volume  of  molecularly 
mixed  fluid  in  a  shear  layer,  provides  a  link  between  the  two  processes. 
Roddam  Narasimha  and  Anatol  Roshko  (private  communications)  conclude, 
somewhat  tentatively,  that  in  the  limit  of  infinite  Reynolds  numbers, 
the  fractional  volume  in  which  dissipation  takes  place,  is  finite.  The 
counterpart  (not  quite  analogously)  in  the  Broadwell-Breidenthal  model 
is  that  in  this  limit,  the  molecularly  mixed  volume  fraction  is  finite 
-  0.3.  Furthermore,  in  the  same  limit,  the  kinetic  energy  dissipation 
rate  is,  in  the  classical  theory,  independent  of  Re  and  the 
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corresponding  assumption  is  made  in  the  model  for  scalar  mixing.  These 
tentative  ideas  will  be  further  pursued  as  they  promise  to  provide  a 
firmer  theoretical  basis  for  the  model  as  well  as  to  lead  to  new  ideas 
about  turbulence  itself. 

With  regard  to  the  role  of  Reynolds  number,  we  have  emphasized  its 
influence  on  mixing  in  gases  (Mungal,  Hermanson  &  Dimotakis  1985)  but  of 
equal  importance,  from  a  fundamental  point  of  view,  is  its  lack  of 
influence  on  mixing  in  liquids.  Koochesfahani  &  Dimotakis  (1986) 
confirm  an  earlier  finding  of  Breidenthal  (1981),  that  there  is  no 
change  in  reaction  product  in  the  Reynolds  number  range  between  about 
2x10^  and  about  8*10^  ( Breidenthal’ s  absolute  measure  of  product  has 
been  shown  to  be  in  error,  but  there  is  no  reason  to  doubt  the  trend,  or 
lack  thereof).  Likewise,  for  high  Re,  turbulent  jet  flame  length  is 
empirically  found  to  be  independent  of  Re,  again  a3  predicted  by  the 
model.  The  practical  inference  to  be  drawn  from  these  analyses  is  that 
the  high  Reynolds  number  asymptotic  range  for  gases  lies  above  -  10^  and 
that,  therefore,  most  (probably  all)  laboratory  combustion  experiments 
are  influenced  by  explicit  molecular  effects.  As  we  have  emphasized 
(Broadwell  &  Breidenthal  1982  and  Broadwell  &  Dimotakis  1986), 
conventional  turbulent  diffusion  theories  do  not  address  these  important 
effects. 

The  axi-symmetrlc  jet  model  is  described  in  the  attached  paper 
(appendix  D),  accepted  for  presentation  at  the  US-France  Joint  Workshop 
on  Turbulent  Reactive  Flows,  6-10  July  1986  (Rouen,  France).  It  is 
suitable  for  engineering  use,  i.e.  the  computations  are  of  practical 
length  with  realistic  chemical  kinetic  systems.  It  will  be  employed  to 
examine  the  "nitric  oxide  problem"  in  methane  and  hydrogen  jets  that  has 
been  discussed  in  earlier  reports. 

With  regard  to  modeling  supersonic  combustion,  there  are  reasons  to 
hope  that  this  model  is  applicable  there  also,  at  least  roughly  so.  The 
Ph.D.  thesis  of  Dimitri  Papamoschou  (1986)  suggests  the  presence  of 
large  scale  structures  in  the  flow,  the  critical  feature  underlying  the 
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model,  and  there  la  no  a  priori  reason  to  doubt  its  applicability.  The 
idea  of  homogeneously  mixed  fluid  is  presently  being  used  as  a  design 
tool  for  the  supersonic  shear  layer;  the  strained  flames  are  not 
explicitly  treated  (see  discussion  is  section  2.4).  We  are  encouraged 
in  this  use  by  the  fact  that  this  form  of  the  model  provided  valuable 
guidance  throughout  the  supersonic  chemical  laser  program,  first 
described  in  Broadwell  (1974).  Other  results  from  that  program  are 
likely  to  be  helpful  in  the  supersonic  combustion  study. 

In  summary,  through  the  work  of  the  past  year,  there  is  promise  of 
showing  that  the  model  is  consistent  with  the  basic  parts  of  classical 
turbulence  theory.  It  has  now  been  put  in  quantitative  form  for  shear 
layers  and  jets,  and  provides  a  unified  interpretation  of  our 
experimental  observations  for  wide  ranges  of  all  the  parameters  that  we 
have  investigated.  J.  E.  Broadwell  is  an  invited  lecturer  on  these 
topics  at  the  NATO  Advanced  Summer  Institute  on:  Disorder  and  Mixing  to 
be  held  at  the  Institut  D'Etudes  Scientif iques  de  Carg^se,  12-27  June 
1987. 


4.0  COMPUTATIONAL  EFFORT 


We  have  tested  several  types  of  Lagrangian  elements  on  a 
three-dimensional  flow  that  features  a  strong  interaction  between  vortex 
structures  —  the  merging  of  two  vortex  rings.  Figure  1 4a  shows  the 
merging  as  computed  by  a  vortex  filament  method  and,  figure  14b,  by  a 
vortex  stick  method.  In  the  filament  method,  tubes  of  vorticity  are 
defined  by  a  sequence  of  material  points  that  move  with  the  local 
velocity.  This  is  an  efficient  representation  for  short  times  but  has 
difficulties  once  the  viscous  reconnection  of  vortex  lines  becomes 
important.  In  the  vortex  stick  method,  vector  elements  of  vorticity  are 
used  that  move  with  the  local  velocity  and  are  stretched  and  rotated  in 
response  to  the  local  strain  field.  This  method  appears  to  have  the 
flexibility  and  efficiency  required  to  handle  vortex  merging  in  three 
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dimensions. 

We  are  continuing  our  efforts  to  reduce  the  computational 
complexity  of  the  vortex  method,  from  N2  to  Nlog(N),  where  N  is  the 
number  of  computational  elements.  The  presently  investigated  scheme, 
proposed  by  Appel,  relies  on  using  a  monopole  (center-of-vorticity ) 
approximation  for  computing  induced  velocities  over  large  distances  and 
a  binary  free  data  structure  to  keep  track  of  which  vortices  are 
sufficiently  clustered  to  warrant  the  clustering  dynamically.  In  a 
separate  investigation  of  mixing  and  chaotic  phenomena,  a  numerical  and 
analytical  study  is  in  progress,  preliminary  results  of  which  were 
recently  reported  at  the  Chaos  87  meeting  in  Monterey  (Leonard, 
Rom-Kedar  4  Wiggins  1987).  Using  non-linear  system  dynamics  methods, 
the  entrainment,  mixing  transition  and  other  dynamical  properties  of  a 
simple  vortex  system  were  analyzed. 


5.0  DIAGNOSTICS,  INSTRUMENTATION  &  DATA  ACQUISITION 


5 . 1  Progress  In  low  light  level  photodetection 

This  part  of  the  work  was  described  as  part  of  section  3-1  • 


5.2  Iwproved  data  acquisition  rate 

We  have  completed  the  first  phase  of  our  efforts  to  improve  the 
maximum  rate  that  we  can  record  data  into  computer  memory  and  or 
computer  disk  to  the  following  numbers. 

a.  Transfers  to  disk  (sustained  for  a  maximum  of  80  MBytes;  at 
1  . 2  MBytes/s. 
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b.  Transfers  to  computer  memory  for  a  maximum  of  4  MBytes  at 
roughly  3.5  MBytes/s. 

c.  Burst  transfers  to  computer  memory  (for  a  maximum  of  1  kByte 
per  burst)  at  5  MBytes/s. 


The  second  part  of  our  effort  to  increase  these  numbers  further  (in 
excess  of  24  MBytes/s)  is  in  progress.  The  design,  based  on  the  high 
speed  VME  bus,  is  nearly  complete.  Fabrication  of  interface  boards  to 
the  new  bus,  designed  for  specific  data  acquisition  applications,  should 
be  in  progress  by  the  end  of  the  summer.  Such  data  rates  should  permit 
real-time  recording  of  two-dimensional  digital  image  data,  as  well  as 
high  framing  rate  linear  array  data  (as  needed  for  the  anticipated 
supersonic  shear  layer  work). 

This  part  of  the  effort  was  co-sponsored  by  the  DoD/URIP  contract 
AF0SR-85-01 53  and  is  performed  by  D.  Lang  in  collaboration  with  P. 
Dimotakis. 
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FIGURE  4.  Low  heat  release,  H2/f2/N0  system  reactor  kinetics.  Solid 
line  for  constant  mass  reactor.  Dot-dashed  line  for 
entrainment-dominated  reactor.  Intercept  of  maximum  slope 
(dashed)  line  with  AT(t)/ATf  -  1  at  t  -  tx-  High  speed 

stream:  8?  H2  ,  0.03?  NO  ,  91.97?  N2  §  *  293  K  & 

Mi  *  O.O63  •  Low  speed  stream:  1  ?  F2  ,  99?  N2  @  T2  *  293  K  & 
M2  ,  0.025  . 
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FIGURE  8. 


Shear  layer  mixing  model  (Dimotakis  1987)  predictions  for 
dependence  of  <5p/5  chemical  product  volume  fraction  on 
Reynolds  number.  Circles  and  squares  from  gas  phase 
measurements  of  Mungal  et  al  (1985).  Triangles  from  liquid 
phase  measurements  of  Koochesfahani  &  Dimotakis  (1986). 
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FIGURE  12.  High  pressure  reactant  vessel  (see  section  3*2) 
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FIGURE  14.  Computation  of  vortex  ring  merging.  (a)  vortex  filaments 
(b)  vortex  sticks. 
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Abstract 

The  linear  spatial  instability 
characteristics  of  Doth  uniform  and  non-uniform 
density  plane  mixing  layers  were  investigated 
talcing  into  account  the  wake  component  of  the 
initial  velocity  profile.  Two  unstable  modes  were 
found.  In  the  shear  layer  mode,  the  growth  of  the 
unstable  disturbance  leads  to  tne  usual 
Kelvln-Helmholtz  roll-up  pattern,  whereas  In  the 
wake  mode,  roll-up  patterns  resemble  those  In  wake 
flows.  It  was  found  that  the  shear  layer  mode 
dominates  the  wake  mode  when  the  density  is 
uniform  across  the  layer.  The  wake  mode,  however, 
can  become  comparable  or  even  stronger  than  the 
shear  layer  mode  If  the  density  of  the  low-speed 
stream  Is  larger  than  that  of  the  high-speed 
stream.  Experimental  evidence  In  support  of  these 
findings  1s  provided. 

Introduction 

The  lnvlscld  linear  instability  of 
two-dimensional  two-stream  plane  mixing  layers  has 
been  studied  extensively  in  the  past.  In  the  case 
of  uniform  density,  Mlchalke1  investigated  the 
single-stream  shear  layer,  while  the  effect  of  the 
velocity  ratio  in  two-stream  mixing  layers  was 
considered  by  Monkewltz  A  Huerre*.  Maslowe  4 
Kelly*  studied  stratified  (non-uniform  density) 
shear  layers  and  showed  that  density  variations 
can  be  destabilizing.  In  all  these  studies,  the 
mean  velocity  profile  has  been  assumed  to  be 
monotonically  increasing  from  the  value  on  the 
low-speed  stream  to  that  on  the  high-speed  stream 
and  usually  the  hyperbolic  tangent  form  is  used. 
It  should  be  noted,  however,  that  under 
experimental  conditions  the  initial  mean  velocity 
profile  almost  always  has  a  wake  component  due  to 
the  boundary  layers  on  the  two  sides  of  the 
splitter  plate.  The  effect  of  the  wake  component 
has  only  recently  come  Into  consideration  with  the 
investigations  of  Mlau'  and  Zhang  et  al.*  for  the 
uniform  density  case. 


•  Present  address:  Department  of  Mechanical 
Engineering,  Michigan  State  University, 

East  Lansing,  Michigan  48824.  Member  AIAA. 

*  Graduate  Student,  Aeronautics. 


The  purpose  of  the  present  worx  :s  to  at „cy 
the  Instability  characteristics  of  both  uniform 
and  non-uniform  density  plane  shear  layers,  wncse 
initial  velocity  profiles  Include  a  wa<e 
component.  The  lnvlscld,  linear,  parallel-flow 
stability  analysis  of  spatially  growing 
disturbances  is  utilized  to  numerically  calculate 
the  range  of  unstable  frequencies  ar.c 
wave-numbers.  The  flow  patterns  resulting  from 
the  amplification  of  the  instability  are  examined 
by  calculating  the  streaklines  and  are  compared 
with  the  experimental  flow  visualization  pictures. 

Formulation  of  the  Problem 

We  consider  the  general  case  of  a  two-stream 
plane  shear  layer  with  U,,  p,  as  tne  free-stream 
velocity  and  density  on  the  high-speed  stream  and 
U,,  p,  as  the  corresponding  quantities  on  the 
low-speed  side  of  the  layer.  All  of  tne 
quantities  used  here  are  normalized  with  tne 
average  velocity  (Ut*U,)/2,  average  density 
(Pi*Pj)/2  and  the  local  layer  thickness  6  as  tne 
length  scale.  Since  we  are  not  aware  of  any  exact 
solutions  for  the  initial  evolution  of  non-uniform 
density  mixing  layers  with  a  wake  component,  we 
assume  that  the  mean  velocity  and  density  profi.es 
have  the  following  forms.  The  mean  velocity 
profile  Is  composed  of  the  usual  hyperbolic 
tangent  profile  plus  a  wake  component  (due  to  the 
splitter  plate)  represented  by  a  Gaussian 
distribution  and  has  the  form 

U(y)  -  1  ♦  lutanh(y)  -  We'y‘ 

where  ly  •  (U,-U,)/(U,»U,)  and  W  is  the  normalized 
wake  deficit.  The  mean  density  profile  has  a 
hyperbolic  tangent  profile  and  is  given  by 

p(y)  -  i  *  lp  tann[(y-y,)/a]  '  l'. 

where  lp  ■  (pi-p»)/'(p,*pl)  and  y,  and  o  adjust  tne 
lateral  position  and  thickness  of  the  density 
profile  relative  to  the  velocity  profile. 

The  disturbance  stream  function  is  written  m 
the  form 

*  -  *(y )  eu«_8t)  •' 

where  a  •  Op^iai  is  the  complex  non-dimensional 
wave-number  and  5  is  the  non-dimensional  frequency 
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which  Is  taken  as  pure  real  for  the  present 
spatial  calculations.  In  the  case  of 

two-dimensional  incompressible  flow  with 
negligible  buoyancy  effects  (i.e.  gravity  is 
Ignored),  it  can  be  shown  that  the  disturbance 
eigenfunction  *  satisfies  the  equation 


♦  *  ♦  (o'/o)  ♦'  *  [«* 


U"  ♦  p  '  U'/p 
U-fl/a 


♦ 


C 


(«) 


where  (  )'  corresponds  to  d/dy.  The  equation 
above  reduces  to  the  Rayleigh  equation  when  the 
density  is  uniform.  A  "shooting”  technique  was 
used  to  solve  this  eigenvalue  equation  with  the 
boundary  conditions 

♦<y  *  *•)  -  e  *ay.  (5) 

Equation  *  was  Integrated  from  both  sides  toward 
y  -  0  and  the  matching  of  e  and  at  this  point 
yielded  the  spatial  growth  rate.  -alt  of  unstable 
disturbances  and  the  corresponding  wave-number, 
<V,  versus  frequency  S. 

The  flow  patterns  resulting  from  the 
amplification  of  the  instability  were  determined 
from  streakline  calculations.  The  procedure  for 
this  calculation,  which  1s  essentially  the  same  as 
that  used  by  Mlcha’ ke‘,  1s  outlined  below.  The 
perturbation  velocities  are  given  by 

u(x,y,t)  -  Real|3e/3y}, 

v(x,y,t)  -  Real{-3S/3x} ,  (6) 

where  *  is  given  by  aquation  3  and  «  is  known  from 
the  solution  of  equation  4  for  a  specific  set  of 
eigenvalues  (a,  3).  The  motion  of  each  fluid 
particle  is  then  given  by 


U 

Fig.  1  Mean  velocity  profiles  for  iifferec 
values  of  H,  lu  .  o. «5  (U,/U,  •  1.38  . 

. W  •  0,  tanh  profile:  -• —  «  • 

-  W  •  0.6:  -  W  •  0.8. 
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dx/dt  «  U(y)  ♦  tu(x.y.t), 

dy/dt  -  e  v(x.y.t),  (7) 

with  the  initial  conditions  x(t-0)  •  x,  and 

y(t-O)  •  yt.  In  equation  7,  U(y)  is  the 

undisturbed  mean  profile  (equation  1)  and  e  is  a 
measure  of  the  initial  magnitude  of  the 

disturbance.  For  the  streaklines  shown  in  the 
present  work,  x(  was  selected  to  be  zero  with 
c  •  0.0005.  Calculations  were  made  by  integrating 
equation  7,  at  different  starting  y  locations, 
forward  in  time  using  the  Euler  method. 


0 


Fig.  2  Instability  characteristics  of  the  mean 
profiles  in  figure  1.  uniform  density.  For  legend 
see  figure  1. 
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Results  and  Discussion 
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The  spatial  instability  character istics  of 
tns  shear  layer  witn  a  wake  ooaponent  were 
calculated  for  a  fixed  velocity  ratio  as  a 
function  of  the  depth  of  tne  wake  deficit  (see  the 
seen  velocity  profiles  in  figure  U.  These 
profiles  can  be  thought  of  as  representing  the 
evolution  of  the  seen  profile  due  to  viscous 
diffusion.  Starting  near  the  splitter  plate  tip, 
a  profile  with  a  large  wake  deficit  develops,  as 
it  convects  downstreaa,  into  or..'  with  no  wake 
component.  The  earn  result,  shown  in  figure  2,  is 
that  when  the  wake  ooaponent  is  present,  there  are 
two  unstable  nodes  as  opposed  to  one  in  tne  case 
of  the  nyperoolic  tangent  profile.  Consistent 
with  previous  results',  as  the  wake  deficit 
increases,  the  neutral  point  of  aode  ?,  the 

stronger  node,  aovea  to  lower  frequencies  and  its 
aaxlaua  amplification  rate  Increases.  We  point 
out  that  the  existence  of  aode  2.  the  weaker  aode, 
has  been  known  fra®  the  work  of  Mlksed*.  His 
results,  however,  were  based  on  teaporal  stability 
calculations.  Note,  froa  figure  2.  that  in  the 
Halt  of  zero  wake  coaponent.  node  <  approaches 
the  tann  profile  solution  while  aode  2  vanishes. 

The  set  of  profiles  in  figure  3  were  used  to 
calculate  tne  behavior  of  the  solution  in  the 

Halt  of  unity  velocity  ratio,  naaely  the  case  of 
pure  wake.  Results  shown  in  figure  *  indicate 
that,  in  tnis  Halt,  tne  two  eodes  of  instability 
still  persist.  In  tne  Halt  of  pure  wake, 

ezaainatlon  of  tne  eigenf  isict  ions  (not  shown  nere) 
revealed  tnat  eodes  '  and  2  approach  tne  ‘sinuous” 
and  ‘varicose*  nodes.  respectively,  of  wake 

instability  (e.g.  see  the  wake  stability 

solutions  of  Mattingly  t  Criminals’). 


Fig.  3  Mean  velocity  profiles 
values  of  iu,  4  m  3.9. 

•  0.0,  pure  wawe: 

-  lu  -  0.30:  -  »u  •  0.4?. 
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Streaklines  were  celculated  in  order  to 
exaalne  the  flow  patterns  that  would  result  froa 
the  amplification  of  the  instable  disturbance  in 
each  of  the  two  nodes  (see  the  previous  section 
for  details).  Calculations  were  perforsed  only 
for  the  case  of  aaxlaua  amplification  rate.  The 
integration  proceeded  until  the  first  structure 
roll-up  appeared.  The  results  are  illustrated  in 
figures  5  and  6.  We  emphasize  that  the  calculated 
streaklines  are  only  intended  to  provide  a 
qualitative  description  of  shear  layer  roll-up 
patterns.  Non-parallel  flow  and  nonlinear  effects 
are  absent  In  these  results. 

The  flow  petterns  in  tne  case  of  tne  tarn 
profile,  figure  5,  are  siallar  to  those  calculated 
by  Mlchalke'  (for  U,-0  shear  layer)  and  snow  tne 
faalllar  roll-up  of  the  snear  layer  into  a  vortex. 
Figure  6  depicts  tne  streaklines  for  tne  two  nodes 
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Fig.  4  Inatabi.lty  cnarecter : st .  ca  .f  -» 
profiles  in  figure  3,  uni  fora  dersit*  r  *■  .  **•- 

see  figure  3- 
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profl.e,  .  3.4?. 


pens  1 ty  tann 


3t-*a«..re*  .»  uniform-density  snear  .ayer 
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if  .ntaoi.ity  -nan  tne  <oa  component  is  present. 
4a  -iota  tnat  tne  amp. if icatlon  of  soda  ’  .eads  to 
tne  uaua.  te./in-Helsnoiti  -  ype  snaar  layer 
-ol.-up.  -nereaa  tne  sods  2  -oil -up  patterns 
-esemo.e  a  -axe  flow.  4a,  tnerefore,  -afar  to 
aodaa  and  2  as  tna  ’snaar  .ayer"  and  "vaxe" 
sodas.  -aspact i  »a,r .  Tha  -axe  soda,  in  tna 

4 r.fors  lansity  :aaa.  is  very  difficult  to  00a erve 
•spar isenta..y  sinca  its  amplification  rata  is 
such  .ass  tnan  tnat  of  tna  snaar  .ayar  soda 
’.nrouffi  scat  of  tna  unstable  frequency  -an (a  saa 
fifui r»  1 
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In  investigating  tna  affact  of  oon-unifors 
lansity.  a  specific  velocity  proflla  -itn 

*  4  •  3  .*8  J 1  ‘J ,  -  3.38  and  a  rasa  deficit  of 
4*38  >as  sa.actad.  4a  a ra  particu.ariy 
intaraatad  . n  tna  caaa  of  .ow-speed  straaa  having 
tna  nifrar  penalty  *he  lansity  proflla  -as 


Fig.  T  Mean  velocity  and  density  pron.es, 
iu  -  0.»5,  l?  •  -0.78  (p,/p,  •  3),  ^0  -  ’6. 

arranfad  to  nava  Its  lnflactlon  point  at  tne 
ilniaua  of  tna  velocity  proflla  and  Its  tlucxness 
auen  smaller  (in  aquation  2,  1/9-16)  tnan  tnat 

of  tna  velocity  proflla,  saa  flgu r«  These 

conditions  ars  sxpactsd  to  nold  in  tne  initia.. 
region  of  the  flow  nesr  tne  splitter  plate  tip. 
The  qualitative  features  of  tna  results  are  not 
sensitive  to  these  conditions  as  long  as  tne 
Penalty  profile  is  "reasonably*  tMn  relative  to 
tne  velocity  profile.  Sea  also  tna  next  section. 
Tabulations  ware  parforsad  for  9,/e,  *  i 
ij  •  -0.78).  “he  ceea  of  nigh-danslty  hign-speea 
straaa  '#,/#,  •  '/8,  i9  .  -o.78)  -as  a.so 

calculated  for  coaparlson. 

Tha  soat  important  resu.  t,  see  figure  S,  is 
tnat  wnsn  tne  nigh  density  is  on  tne  low-speea 
side,  tne  two  instability  modes  nave  sisi.ar 
aapllf icatlon  rates.  In  fact,  tna  normally  -eax 
-ana  soda  in  tna  case  of  unifora  density,  now  ras 
a  silently  larger  growtn  rata  tnan  tne  snear  layer 
mode.  The  corresponding  streamlines  at  maximum 
aapl if Icatlon  rate)  for  these  two  sodas,  figure  ?, 
again  illustrate  the  shear  layer  and  -axe  type 
-oli-up  patterna.  4a  also  note,  from  figure  8, 
tnat  naving  tna  nigh  density  on  tna  nign-speec 
straaa  Poes  not  altar  the  relative  significance  pf 
tna  two  instability  sodas  compared  to  the 
uniform-density  shear  layer. 

Tha  findings  described  above  imply  tnat. 
depending  on  tne  speccrua  of  tna  disturbances 
tna  flow  and  tna  extant  of  tna  persistence  of  tne 
-axe  coaponent  in  the  downstreaa  region,  a  snear 
layer  of  non-unifora  penalty  say  not  -oil  up  lima 
tna  usual  Kalvin-Helaholtz  structures  Out  more 
lixe  a  wane.  These  results  also  suggest  tnat, 
under  tna  right  flow  conditions,  both  snaar  layer 
and  -axe  sodas  of  instability  may  exist 

s i aul taneoua. y  and  interact  witn  «acn  other. 


Sine*  the  shear  layer  and  wake  modes  can  be 
equally  strong  when  the  low-speed  stream  has  the 
high-denslty  fluid,  an  attempt  was  made  to 
experimentally  observe  these  modes.  A  shear  layer 
(U,/U,  •  0.38)  between  a  high-speed  stream  of 
Helium  and  low-speed  stream  of  Argon  (p,/pi  •  10) 
was  forced  acoustically.  The  flow  visualization 
by  Schlieren  photography,  figure  10,  shows  that 
both  shear  layer  and  wake  modes  can  be  generated 
in  a  two-stream  mixing  layer.  Note  that  the  wake 
mode  pattern,  in  figure  10(b),  appears  to  approach 
that  of  the  shear  layer  mode  toward  the  right  side 
of  the  photograph.  This  is  to  be  expected,  since 
as  the  flow  moves  downstream  the  wake  component  of 
the  velocity  profile  ultimately  vanishes  and  only 
the  shear  layer  instability  mode  remains. 


Effect  of  the  Density  Profile  Thickness 


The  case  of  uniform  density  can  be  considered 
to  be  equivalent  to  that  of  non-uniform  density 
with  a  very  large  (in  fact,  infinite)  thickness 
relative  to  the  mean  velocity  profile.  A 
question,  then,  remains  as  to  how  a  weak  wake  mode 
in  a  uniform-density  shear  layer  transforms  into 
one  equally  strong  as  the  shear  layer  mode  when 
the  high-denslty  fluid  is  carried  on  the  low-speed 
side.  To  shed  some  light  on  this  question,  the 
instability  characteristics  of  the  non-uniform 
density  shear  layer  were  calculated  as  a  function 
of  the  density  interface  thickness  (l.e.  a  was 
varied  in  equation  2).  The  mean  velocity  and 
sample  density  profiles  are  shown  in  figure  11. 
The  results,  see  figure  12,  illustrate  that  the 
density  profile  thickness  must  be  smaller  than  a 


Fig.  11  Mean  velocity  and 
different  values  of  a, 
Ip  •  -0.78  (pa/p,  •  8). 
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certain  value  before  the  wake  mode  becomes 
dominant.  For  these  particular  profiles,  for 
example,  it  is  required  that  the  thickness  of  the 
density  profile  be  at  least  2.9  times  smaller  than 
the  thickness  of  the  velocity  profile.  On  the 
other  hand,  when  the  high-denslty  fluid  is  on  the 
high-speed  side,  regardless  of  the  density  profile 
thickness  (lowest  value  calculated  was  a  •  1/16), 
the  shear  layer  mode  of  instability  was  found  to 
be  always  dominant. 


An  Interesting  feature  in  figure  12  is  that 
the  appearance  of  a  strong  wake  mode,  between 
a  •  1/2.8  and  1/2.9,  seems  to  be  a  resonance 
phenomenon.  This  "repelling"  of  two  otherwise 
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Fig.  13  Shape  of  ( U"  *  p'U’/p)  profile  for 

different  denalty  profile  thicknesses. 

•••••  p  •  1;  -•-•-  q  m  1/2;  — —  9  *  1/2.5; 

—  9  *  1 /3. 

Identical  eigenvalues  la  tnought  to  t>«  alallar  to 
the  "Eckart"  raaonanca*.  which  la  the 
quant urn- mechanical  analog  of  the  problem  of  two 
potential  minima  and  the  penetration  of  the 
potential  barrier  between  thee.  In  coaparlaon 
with  the  Eckart  reaonan ce  phenomenon.  It  aeena 
that  the  behavior  of  the  tare  (W  ♦  p'U'/p)  in 
equation  a,  in  particular  the  preaence  of  two 
maxima  in  lta  profile,  la  the  determining  factor. 
A  plot  of  thla  term  aa  a  function  of  o.  aee 
figure  13,  ahowa  the  appearance  of  two  maxima 
cloae  to  the  value  of  o  where  the  reaonance  occur a 
(aee  figure  12) . 

Conclusions 

The  inatablllty  character l at lea  of  uniform 
and  non-uniform  denalty  plane  ahear  layera  were 
Inveatigated.  The  mean  velocity  profile  included 
a  wake  component  In  order  to  take  account  of  the 
effect  of  the  boundary  layera  on  the  two  aldea  of 
the  aplltter  plate.  The  range  of  unatable 
frequenclea  and  wave-numbera  and  the  flow  patterna 
reaultlng  from  the  amplification  of  the 
Inatablllty  were  calculated  ualng  the  inviacld, 
linear,  parallel-flow,  apatlal  atablllty  analyala. 

It  wax  found  that  the  ahear  layer  with  a  wake 
component  haa  two  unatable  modea.  The  growth  of 
tne  unatable  diaturbance.  In  the  ahear  layer  mode, 
reeulte  In  the  uaual  Kelvin-Helmholtx  roll-up 
patterna,  while  In  the  wake  mode,  the  flow 
patterna  reaemble  a  wake  atructure.  When  the 
denalty  la  uniform,  the  amplification  rate  of  the 
wake  mode  la  dominated  by  that  of  the  ahear  layer 
mode.  If  the  low-apeed  stream  carries  the 
hlgh-denslty  fluid,  however,  the  two  modea  can 
become  comparable  In  amplification  rate.  The 
oppoelte  arrangement  with  the  hlgh-denslty  fluid 
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on  tne  hign-apeed  side  denaves  sisi.ar  to  t.-.e 
uniform  density  case  m  t.nat  tne  shear  .ayer  mode 
remains  the  dominant  node.  3c.nl  :eren  flew 
visualization  pictures  of  a  snear  layer  petween  i 
high-speed  stream  of  a  light  gas  Helium,1  ar.d  a 
low-speed  stream  of  a  neavy  gas  (Argon,  confirmed 
tnat  both  modes  of  instability  could  se  excited 
experimentally.  For  the  waxe  mode  to  become 
dominant,  tne  thickness  of  the  density  profile, 
relative  to  the  velocity  profile  thickness,  must 
be  smaller  than  a  certain  value.  The  onset  of  a 
strong  wake  mode,  aa  tne  density  profile  thickness 
Is  reduced,  appears  to  exhibit  a  resonance 
phenomenon.  It  la  believed  that  this  oenav-.cr  is 
similar  to  the  "Eckart"  resonance. 


Aotcaowledmsmsnta 


We  are  greatly  Indebted  to  Professor  Tosr.i 
Kubota  who  is  an  unofficial  co-autnor  of  tnis 
paper.  He  always  had  the  answers  to  our  questions 
and  the  time  to  discuss  Ideas.  Alio,  tne  constant 
interest  In  thla  problem  by  Professor  Paul 
Dlmotakla  and  the  late  Professor  Lester  Lees  is 
much  appreciated.  This  work  was  supported  by  tne 
Air  Force  Office  of  Scientific  Reaearcn  Grant  Ho. 
AF0SR-83-021 3. 


Rer 


1.  Michalke,  A.,  "On  Spatially  Growing 
Disturbances  In  an  Inviacld  Shear  Layer.", 
fluid  Mach.  23(3).  i 965.  pp.  521  -544. 


2.  Monkewltz,  P.  A.  and  Huerre,  P. ,  "Influence 
of  the  Velocity  Ratio  on  tne  Spatial 
Instability  of  Nixing  Layers.",  Phys.  Fluids 
25(7).  1982,  pp.  1137-11*3. 


3.  Haalowe,  S.  A.  and  Kelly,  R.  E. ,  "Inviscid 
Inatablllty  of  an  Unbounded  Heterogeneous 
Shear  Layer.",  Fluid  Nech.  *8(2: ,  '9"', 

pp.  405-«l  5. 


4.  Nlau,  J-J.,  "An  Experimental  Study  on  tne 
Instability  of  a  Mixing-Layer  vitn  Laminar 
Wake  aa  the  Initial  Condition.",  Ph.G. 
thesis.  Brown  University,  '984. 


Zhang,  Y-Q,  Ho,  C-H.  and  Monkewitz,  P.,  "The 
Mixing  Layer  Forced  by  Fundamental  and 
Subharmonlc.",  LAMINAR-TURBULENT  TRANSITION. 
Proc.  IUTAH  Symp . ,  Novosibirsk,  USSR.  July 
9-13.  1984.  Sprlnger-Ver lag,  pp.  385*395. 


6.  Mlksad,  R.  W.,  "Experiments  on  the  Nonlinear 
Stages  of  Free-Shear-Layer  Transition.",  J^ 
Fluid  Mech.  56(4).  1972,  pp.  895-T'9. 


7.  Mattingly,  0.  E.  and  Criminals,  W.  0..  "The 
Stability  of  an  Incompressible  Two-Dimensional 
Wake.",  Fluid  Mech,  51  (2) .  '972.  pp . 

233-272. 


8.  Eckart,  C.,  "Internal  Waves  in  the  Ocean.". 
Phya.  Fluids,  4(7).  1961,  pp.  79i-?99. 


m-v -warm 


& 


VJ 


".V 


.*• 


£ 


VJ 


The  following  represents  a  listing  of  the  Caltech  Pressure  Vessel 
Specifications : 

1.  The  internal  tank  volume  must  be  in  the  range  of  1.2  m3  to 
1.3  m3  (42  ft3  to  44  ft3). 

2.  The  working  pressure  range  is  vacuum  to  10^  Pa  (vacuum  to 
1500  psi) . 

3-  The  working  temperature  range  is  273  K  to  615  K  (32°F  to 
650°F) . 

4.  The  working  gases  are  (dry)  H2,  NO,  N2.  He,  Ar. 

5.  The  tank  material  is  carbon  steel.  The  preferred  type  is 
SA5 16-70. 

6.  The  tank  geometry  is  cylindrical  with  2:1  semi-elliptical  end 
caps  (refer  to  the  attached  sketch  for  clarification). 

7.  One  of  the  end  caps  shall  be  removable.  This  will  require  a 
flanged  connection  as  depicted  in  the  attached  sketch. 

8.  Nominal  access  to  the  tank  is  via  two  0.1  m  (4.0  in)  diameter 
weld  neck  flanges,  one  at  either  end. 

9.  The  outer  diameter  of  the  tank  is  0.91  ®  (3-0  ft). 

10.  The  maximum  flange  face  to  flange  face  height  of  the  tank  is 
2.75  ®  (9-0  ft). 

11.  The  interior  of  the  tank  shall  be  packed  in  the  following 

manner:  Vertical  rolls  of  aluminum  mesh  screen  shall  fill  the 

cylindrical  portion  of  the  tank,  except  for  the  cross-section 
occupied  by  a  central  alignement  tube  and  necessary  clearances 
between  the  packing  rolls  and  the  internal  tank  walls. 

12.  The  aluminum  rolls  are  to  be  supported  by  three  stainless  steel 
grates,  one  at  either  end  and  one  in  the  middle.  These  grates 
are  to  be  bolted  to  mounting  rings  which  are  in  turn  welded  to 
the  inner  wall  of  the  tank. 

13.  Caltech  will  supply  the  aluminum  screen  and  the  central 
alignment  tube.  The  tank  manufacturer  will  be  responsible  for 
the  packing  supports  and  the  overall  assembly. 

14.  The  tank  is  to  be  mounted  vertically  on  top  of  a  skirt.  The 
skirt  must  have  four  0.45  m  (18  in)  circular  holes  cut  into  it 
aligned  in  a  cross  fashion.  The  skirt  will  have  a  connecting 
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ring  welded  to  its  bottom  which  is  suitable  for  bolting  the 
tank  assembly  to  the  floor.  The  height  of  the  centers  of  these 
holes  is  to  be  agreed  upon  at  the  time  of  manufacture. 

15.  The  exterior  of  the  tank  is  to  be  wrapped  with  12  electrical 
heating  tape  elements.  Separate  electrical  connector  pairs 
shall  be  provided  for  each  element  through  the  insulation  (see 
below) .  Two  thermocouples  will  be  mounted  for  each  heating 
element  (for  a  total  of  24)  on  the  tank  exterior  wall,  with 
electrical  connection  pairs  again  made  accessible  through  the 
insulation  material.  Caltech  will  provide  the  heating  tape 
elements  and  thermocouples,  which  will  be  selected  in 
consulation  with  the  tank  manufacturer.  Total  heating  power 
shall  be  in  the  range  8  kW  to  10  kW. 

16.  The  exterior  of  the  tank  is  to  be  insulated  in  such  a  fashion 
that  there  are  no  exposed  surfaces.  The  outside  temperature  of 
the  insulation  shall  not  exceed  ambient  room  temperature  by 
more  than  10  K  (18°F)  when  the  tank  interior  is  at  its  maximum 
temperature  of  615  K  (650°F). 

17.  The  tank  assembly  will  be  stress  relieved. 

18.  The  interior  of  the  tank  is  to  be  sandblasted  clean  after 
stress  relieving. 

19.  The  hydro  test  is  to  be  conducted  in  the  prescence  of  Caltech 

representatives  and  will  follow  a  mutually  agreed  upon 

procedure , 

20.  After  final  assembly,  the  tank  is  to  be  purged  with  nitrogen 
and  shipped  with  a  slight  nitrogen  atmosphere  overpressure. 

21.  The  tank  manufacturer  shall  be  responsible  for  shipping  costs 
to  Caltech. 

22.  The  tank  is  to  be  certified  by  the  manufacturer  according  to 
the  ASME  code  and  any  applicable  local  codes.  In  particular, 
this  certification  must  cover: 

a.  tank  operation  at  the  specified  temperature  and  pressure 
range ,  and 

b.  the  use  of  hydrogen  as  a  working  gas  at  the  specified 
maximum  temperature. 

In  addition, 

c.  the  tank,  the  skirt  support  structure  and  the  mounting 
provisions  shall  be  certified  by  the  manufacturer  to  one 
and  a  half  times  the  local  earthquake  code  specifications. 
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Laser  Soot-Scattering  Imaging  of  a  Large  Buoyant  Diffusion  Flame 
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A  novel  diagnostic  technique,  which  makes  use  of  laser  light  scattered  by  soot  particles,  was  used  in  an  effort  to 
identify  the  flame  sheets  within  a  natural  gas  diffusion  flame.  Soot  particles,  inherently  created  and  consumed  in  the 
flame,  were  used  as  the  scattering  medium,  which  obviated  the  need  for  externally  supplied  seed  material.  Since  no 
foreign  material  was  added  to  the  flame,  the  current  technique  can  be  considered  truly  nomntrusive.  The  soot 
distribution  within  a  large  buoyant  natural  gas  diffusion  flame  is  argued  to  be  a  reasonable  marker  of  the  flame  sheets 
Measurements  made  in  47  4-190  kW  natural  gas  flames  stabilized  on  a  0.5  m  diameter  burner  show  that  the  flame 
sheets  are  highly  wrinkled  and  convoluted  surfaces.  The  flame  sheets  are  distributed  fairly  uniformly  within  the 
instantaneous  volume  of  the  flame,  based  on  images  of  the  associated  soot,  and  the  instantaneous  tlame  volume  is 
devoid  of  soot  for  40-60%  of  the  time.  When  soot  is  present,  it  is  observed  as  thin  sheets  which  become  narrower  in 
regions  where  the  average  strain  rate  is  estimated  to  be  greater. 


INTRODUCTION 

The  complex  internal  flow  field  produced  by  a 
large  diffusion  flame  in  the  region  in  which  rapid 
heat  release  occurs  is  not  yet  well  understood. 
Visual  studies  of  the  flame  luminosity  have  been 
useful  in  studying  some  phenomena,  such  as  the 
flame  height  and  the  unsteady  ‘‘puffing-’  charac¬ 
terizing  many  of  these  flames  [1,  2).  However, 
these  methods  do  not  allow  visualization  of  the 
often  complex  internal  flows  which  are  relevant  to 
the  entrainment  rate  and  heat  release  rate  in  the 
flame,  nor  do  they  suggest  a  sufficiently  detailed 
model  required  for  calculation  of  radiant  flux  to 
the  fuel  beds.  More  recently,  laser  lighting  has 
been  used  effectively  to  study  flows  in  which 
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combustion  is  an  important  feature  [3].  In  these 
cases,  particles  are  used  as  markers  for  the  fuel, 
oxidizer,  or  products.  The  particles  may  be 
created  or  destroyed  in  the  reaction,  or  may  be 
unaffected  by  the  combustion  region.  These  meth¬ 
ods  require  either  upstream  seeding  or  auxiliary 
chemical  reactions  to  produce  the  particles. 

The  purpose  of  the  current  investigation  was  the 
demonstration  of  the  laser  soot-scattering  tech¬ 
nique  as  a  useful  tool  in  studying  the  structure  of 
the  flame  and  the  location  and  nature  of  the 
instantaneous  combustion  interfaces.  Since  soot 
that  is  created  and  consumed  in  the  flame  itself  is 
used  as  the  scattering  medium,  no  extraneous 
seeding  or  chemistry  is  required.  The  technique 
has  been  applied  in  a  study  of  large  buoyant 
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diffusion  flames  with  the  intent  of  elucidating  the 
complex  internal  structure  of  the  large  vortical 
pufflike  features  that  dominate  the  heat  release 
region  of  these  flames. 

QUALITATIVE  DESCRIPTION  OF  THE 
FLAME 

The  flame  studied  in  this  investigation  models 
liquid  pool  fires  and  related  buoyant  flames  by  a 
bed  of  glass  beads  with  fuel  supplied  by  the 
injection  of  natural  gas.  The  associated  flow  is 
highly  unsteady  and  buoyancy-dominated,  and  the 
heat  release  regions  of  the  resulting  flames  are 
made  visible  by  soot  formed  in  the  combustion 
processes.  The  transition  between  laminar  and 
turbulent  flows  and  the  operational  distinction 
between  unsteady  laminar  flows,  which  contain 
highly  wrinkled  laminar  flames,  and  turbulent 
flows  is  often  obscure:  a  body  of  data  does  not  yet 
exist  from  which  a  convincing  set  of  criteria  can  be 
developed  to  define  the  lammar-to-turbulent  tran¬ 
sition  in  these  flames. 

For  flames  with  a  height-to-source  diameter 
ratio  (ZfID)  greater  than  1.0  and  for  which  the 
initial  momentum  flux  is  negligible  compared  with 
the  buoyancy  force,  the  flame  surface  pinches  in 
periodically  toward  the  axis  of  symmetry  just 
above  the  origin  of  the  flame  and  produces  a 
distinct  structure  in  the  flame.  This  structure  has 
the  appearance  of  a  large  irregular  donut-shaped 
vortex  ring,  which  rises  slowly  above  the  source 
and  defines  the  top  of  the  flame  when  the  fuel  it 
contains  burns  out.  In  a  flame  with  a  large  height- 
to-diameter  ratio,  several  of  these  structures  are 
present  at  any  given  time,  and  for  flames  with  a 
height-to-dtameter  ratio  of  two  to  three,  only  one 
structure  is  visible  at  a  time.  For  flames  with 
height-to-diameter  ratio  less  than  one,  the  process 
cannot  be  seen  in  the  visible  flame,  but  does 
appear  in  the  shadowgraph  images  of  the  rising 
plume  of  hot  products.  For  the  flames  in  the 
present  study,  I  <  Z,lD  <  3. 

The  puffing  process  has  been  the  subject  of 
numerous  investigations,  e  g.,  Thomas  et  al.  (I), 
and  recent  reviews  are  given  in  Zukoski  et  al  (4). 
Beyler  [5],  and  Cetegen  et  al.  (6)  The  dominant 
frequency  appears  to  scale  as  \g/D,  for  a  wide 


range  of  diameters  and  fuels.  Here  g  is  the 
gravitational  constant  and  D  is  the  burner  diame¬ 
ter  The  fluctuations  are  most  marked  when  the 
flame  height-to-diameter  ratio.  Z.  D,  is  between 
one  and  three,  and  the  amplitude  of  the  fluctua¬ 
tions  decreases  markedly  as  Z,<  D  increases  from  3 
to  20. 

The  large  size  of  the  structures  and  the  speed  at 
which  they  evolve  suggest  that  they  play  an 
important  part  in  fixing  the  entrainment  and  heat 
release  in  the  lower  part  of  the  flame. 

There  have  also  been  a  number  of  studies  of  the 
radiation  flux  produced  by  buoyancy-controlled 
diffusion  flames.  With  the  exception  of  a  tech¬ 
nique  developed  by  Markstein  [7],  these  experi¬ 
ments  are  based  on  absorption  and  emission 
measurements  of  the  flame  made  along  a  line  of 
sight  through  the  flame,  and  thus  do  not  give  any 
information  concerning  the  structure  of  the  inte¬ 
rior  of  the  flame.  Markstein  s  technique  uses  a 
probe  which  had  a  resolution  of  about  60  mm.  In 
contrast,  the  present  technique  is  nonintrusive  and 
can  be  used  to  give  qualitative  data  for  a  cubical 
resolution  element  with  a  dimension  of  1-2  mm  on 
a  side.  However,  it  will  not  give  the  quantitative 
data  required  to  determine  flame  radiation  that 
Markstein's  intrusive  probe  will  give. 

DESCRIPTION  OF  APPARATUS 

A  50  cm  diameter  burner  was  used  to  deliver 
natural  gas  fuel  to  the  flame.  The  fuel  passed 
through  a  4  cm  deep  porous  bed  of  2  mm  spherical 
glass  beads  whose  surface  was  made  flat  and  flush 
with  the  outer  metal  edge  of  the  burner  The  sides 
of  the  burner  were  vertical  and  extended  to  the 
floor  located  an  adjustable  distance  from  the 
burner  surface,  but  always  maintained  at  greater 
than  one  burner  diameter  distance  in  order  to 
minimize  ground  interference. 

The  burner  was  located  in  the  center  of  a  2.4  m- 
area  enclosed  on  the  sides  by  a  double  layer  of 
aluminum  screen  used  to  reduce  the  strength  ot 
flow  disturbances  present  in  the  laboratory  air 
The  screens  were  2  4m  high  and  above  them  was 
a  2  4  m;  hood  which  removed  the  exhaust  gases 
from  the  laboratory. 

A  schematic  of  the  optical  arrangement  is  shown 
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LASER  SOOT-SCATTERING  IMAGE 

in  Fig.  I.A3WAr"  laser  beam  <514.5  nmt  was 
passed  through  a  steering  and  focusing  apparatus 
which  allowed  accurate  alignment  to  the  axis  of 
the  line  segment  in  the  flame  that  was  to  be  imaged 
and  permitted  the  imaged  segment  to  be  translated 
vertically  with  a  high  degree  of  parallelism.  A 
spherical  lens  collimated  the  beam  to  a  long  thin 
waist  within  the  flame  and  a  polarization  rotator 
maximized  the  scattered  intensity  in  the  viewing 
direction.  The  r-direction  is  taken  along  the  line 
from  the  periscope  to  the  beamstop  and  the  z- 
direction  is  taken  vertically  from  the  face  of  the 
burner,  with  the  center  of  the  burner  face  defined 
as  the  origin. 

The  collection  optics  for  the  scattered  light 
consisted  of  a  heat  reflecting  filter,  one  of  several 
narrow  band  interference  filters,  an  imaging  lens, 
and  a  linear  solid-state  camera.  This  camera  used  a 
(RETICON)  1024  element  photodiode  array  to 
image  the  illuminated  line  segment  up  to  586  times 
per  second.  The  interference  filter  was  used  to 
pass  only  that  light  within  a  narrow  band  <  1- 10  nm 
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F 1 I  The  -  i  mm  Jut  meter  beam  from  j  *  W  Ar  aser 
1 5 14  5  nm  Imei  was  Jirected  ihrtKjgh  ihc  buosjnt  diffusion 
llame  above  the  50  cm  burner  A  segment  of  me  ►vram  was 
imaged  on  a  linear  photodiode  arrav  1 1024  elements)  docked  at 
300  or  600  kHz  Simultaneous  high-speed  ill’  tpsi  movies 
incorporating  two  orthogonal  views  m  the  lumintHis  soot  could 
be  taken  in  synchrony  with  the  linear  arrav  data  The  c 
direction  was  vertical  and  the  r  direction  was  parallel  to  the 
burner  diameter 


at  514.5  nm)  centered  at  the  wavelength  of  the 
incident  laser  light.  Given  the  narrow  bandwidths 
of  these  filters,  very  little  luminosity  was  passed 
through  to  the  camera. 

Because  the  flow  evolved  significantly  as  it 
passed  the  measurement  station  lie.,  relative 
velocities  within  a  structure  were  comparable  to 
the  celerity  of  the  structure),  the  distinction 
between  an  r-t  diagram  and  dn  r-z  "snapshot" 
was  especially  important  1  Structures  were  in  the 
process  of  changing  dramatically  as  they  crossed 
the  measured  line  segment  so  that  the  signature 
they  left  on  an  r-t  diagram  was  quite  distorted 
relative  to  their  appearance  at  an  instant  in  time 
To  make  a  correspondence  between  the  r-t  data 
and  an  instantaneous  view,  a  high-speed  movie 
camera  was  synchronized  to  photograph  two 
perpendicular  views  of  the  flame  luminosity  simul¬ 
taneously  with  the  photodiode  array  images  ill' 
frames  per  second  or  one  movie  trame  for  every 
2.5  or  5  photodiode  scans  i  The  scattered  laser 
radiation  could  not  be  observed  consistently  in  the 
high-speed  movie  films,  although  the  imaged  axes 
could  be  discerned 

A  computerized  data  acquisition  system  based 
on  a  DEC  PDP  1  i  23  CPU  was  used  to  coordinate 
data  collection  The  photodiode  array  was  inter¬ 
faced  to  this  system  through  a  high-speed  4  D 
convener  designed  and  built  at  Caltech  This 
computer  data  acquisition  system  has  been  de¬ 
scribed  previously  |8| 

Both  an  r-i  diagram  and  a  snapshot  an  be  viewed  as  a 
two  dimensional  slice  through  an  ;nherentlv  tour-dimen- 
Monal  obieci  (three  spatial  dimensions  plus  time'  A  snapshot 
saptures  two  spatial  dimensions,  while  an  r - (  diagram 
v aptures  one  spatial  dimension  and  time  The  two  representa¬ 
tions  offer  different  i  redact  ionaf  i  views  « *  f  the  entire  'Nest 
and  >ni\  ,n  the  context  of  a  given  set  >?  questions  can  either 
considered  superior  V  movie  *r  j  two-dimensional 
electronic  imaging  time  senes  offers  j  ihree  dimensional 
own  spatial  plus  time)  view  and  clearly  contains  more 
information  than  either  two-dimensional  view  if  the  spatial 
and  temporal  resolutions  are  comparable  in  all  cases  In  the 
present  investigation  neither  a  snapshot  nor  anv  other  two- 
dimensional  imaging  of  the  soot  scattering  was  possible  with 
ihe  present  illumination  source  i3  W  Ar*  laser),  given  the 
high  background  due  to  the  flame  luminosity  and  due  to  the 
drop  m  intensity  in  forming  a  sheet  ol  light  A  much  higher 
power  laser  or  a  pulsed  laser  and  a  last -shuttered  detector 
would  be  required  to  overcome  this  experimental  difficulty 


12 


RICHARD  C.  MIAKE-LYE  and  STEPHEN  J.  TONER 


RESULTS  AND  DISCUSSION 

Data  were  collected  at  various  z/Z>  values,  from 
0.1  to  2  0.  and  for  various  values  of  the  heat 
releases,  from  47.4  to  190  kW,  as  shown  in  Table 
1  Inlet  Richardson  numbers  for  the  fuel  based  on 
the  burner  diameter.  ( t\pgD)/(p(J :),  varied  from 
about  100.000  at  the  lower  heat  rates  to  5000  at  the 
higher  rates.  The  Reynolds  numbers  at  the  inlet 
ranged  from  100  to  600.  but  increased  down¬ 
stream  as  the  velocities  increased  due  to  buoyancy. 
For  each  entry  in  Table  1  marked  with  x  ,  four 
sets  of  r-t  diagram  data  wre  collected  at  a  600  kHz 
docking  rate.  One  of  these  four  data  sets  was 
collected  in  synchrony  with  the  high-speed  movie 
for  each  entry.  For  the  entries  marked  with  an 

asterisk,  an  additional  data  set  was  also  collected 
at  a  300  kHz  clocking  rate  synchronized  with  the 

table  i 

Flame  Lengths  and  Axial  Measurement  Stations'1 
<*>  47  4  67  0  95  0  134  190  kW 

Flame  Lengths 

/  D  13  16  I."?  2,3  2.82 

A  xml  Position 

:  D 

»>  . 

>1  J  X 

I)  ^ 

I)  * 

0  -5 

i  i) 

:  o 

'  The  experimental  conditions  that  were  studied  are  indicated  in 
tne  table  The  columns  are  headed  by  the  heat  release  rate  <or 
fuel  flow  rale)  in  kW  The  rows  are  labeled  with  the  z  D 
values  For  ait  the  cases  marked  with  an  *  four  sets  or 
0  6?  s  Juration  r-t  data  at  600  kHz  clocking  rate  were 
collected  One  set  was  taken  m  synchrony  with  a  high  speed 
movie  For  those  entries  marked  with  an  asterisk  A  Hot  '00 
kHz  clocking  rate  data  were  also  taken  m  synchrony  with  a 
high  speed  movie  The  M)0  kHz  data  sets  were  image 
prrxessed  and.  for  the  entries  in  the  1^0  kW  column  digital 
mages  are  shown  in  Figs  2-8  The  nondimensional  f.ame 
length.  Z.  D  based  on  50^  miermittency  of  the  tlame.  is 
given  for  each  heat  release  rate  studied 


high-speed  movie.  The  slower  clocking  rate  al¬ 
lowed  nearly  three  times  more  data  to  be  collected 
due  to  computer  data  acquisition  limitations,  albeit 
at  a  lower  temporal  resolution 

For  the  600  kHz  data,  core  memory  limited  the 
data  collection  to  384  scans  of  512  pixels.  This 
amounted  to  0.66  s  of  data  collection.  All  1024 
pixels  were  read  out:  only  the  central  512  were 
recorded  for  ;/ D  <  0.5,  while  every  alternate 
pixel  (again  512  total)  was  recorded  for  ;/D  ^ 
0.5.  For  the  300  kHz  data,  core  memory  was  not  a 
limitation  since  the  data  could  be  written  directly 
to  disk  at  this  speed,  so  that  1024  scans  of  512 
pixels  (or  more)  could  be  accommodated.  Thus 
these  data  sets  represent  3  5  s  of  data  collection 
Since  the  fire  puffs  at  about  2  Hz.  the  higher 
clocking  rate  data  correspond  to  about  one  puffing 
cycle,  while  the  slower  clocking  rate  data  corres¬ 
pond  to  about  seven  cycles 

Flame  lengths  (61.  Z,iD ,  for  various  heat 
releases  studied  here  are  also  shown  in  Table  I 
These  flame  lengths  are  based  on  50%  intermit- 
tency  ot  the  tlame  The  data  presented  in  this 
paper  correspond  to  the  highest  heat  release  values 
collected  with  a  300  kHz  clocking  rate  Except  tor 
the  major  difference  of  tlame  length  (Z.  D)  as  seen 
in  Table  1 .  the  results  presented  are  representative 
of  all  heat  release  values  where  Z.  D  consistently 
reached  the  given  axial  stafon.  i  e  :  <  Z. 

Figures  2-8  show  r  t  diagrams  tor  a  heal 
release  of  190  kW  and  for  c  D  values  ot  <)  1 .  n  2. 
0  35.  0  5.  0  '5.  I  0,  and  2  0  taken  with  j  300  kH/ 
clocking  rale  Darker  regions  indicate  larger 
values  of  scattered  intensity,  although  no  effort  ha> 
been  made  to  interpret  this  intensity  quantitative^ 
since  both  the  particle  size  distribution  and  their 
number  density  affect  ihe  scattered  intensitv 
in  Fig  2.  vertical  lines  indicate  the  locations  .it 
the  pixels  corresponding  to  ihe  center  ot  the  burner 
and  the  burner  edge  'These  lines  were  added  to 
the  image  digitally  after  processing  ihe  Jala  1  The 
same  optical  geometrv  was  used  tor  Figs  2  4  In 
Figs  5-8  a  different  optical  geometrv  thai  itnaced 
more  ihan  a  full  diameter  ot  the  burner  wjs  used 
The  two  outer  edges  of  the  burner  are  marked  ii-r 
this  geometry  in  Fig  5  by  vertical  lines  added  n 
the  same  manner  as  those  in  Fig  2 

In  Fig  3  tour  sections  through  the  imjge  are 
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memuttenc v  which  allows  4  mine  quanuiadve 
anaivso  it  hesc  ideas  and  will  tie  discussed 
bell  iw 

A'  previousiv  stated  the  high  speed  movies 
liter  j  vompirmentarv  view  >1  'he  flan*  struc 
’ure  The  movie  camera  images  the  luminous  soot 
iver  the  whole  flame  not  lust  (he  sout  intersecting 
■ne  aser  beam  Also  the  movie  captures  most  the 
flame  region  at  an  instant  m  tune  >0  that  the  flame 
structure  an  hr  seen  ilirectlv  . n  ir  ;)  space  A 
sequence  >t  high  speed  movie  half  frames  tot  the 
:  (1  -  '  1  *•  A#)  k\A  h*)  kHz  Jala  set  s  sho  1  ;n 
Fig  .11  I  he  .  orresponding  r  1  diagram  is  Fig  *• 

The  putting  behavior  < it  the  flame  s  evident  m 
both  'he  '  1  diagrams  jnd  -n  the  high  speed 
movies  Sear  'he  burner  a  strong  H/  periodic  its 
s  evident  n  the  r  1  data  although  individual 
veparated  putts  have  riot  formed  at  this  stage  Ai 

U  <t  "  J*  and  he’,  ond  Jistinc'  put*-  an  he 

ihscrvcd  Hv  ntegrating  at  toss  '  and  plotting  ihc 


tesulling  average  intensitv  as  a  function  ot  lime 
the  .*  Hi  periodic  il\  an  he  easilv  displaved 
graphivallv  as  is  illustrated  in  Fig  1  I  This  r  pe  it 
phx  ilearlv  shows  that  the  putting  is  a  c|uasi 
periodu  phenomenon  and  that  the  period  ot  the 
vwlit  behavior  is  variable  within  vome  range  The 
mean  frequence  t.  of  the  puffing  Jepends  little  it 
at  all  on  the  heal  release  rate  and  has  a  nondimen 
sionali/ed  value  of  /  ,D  g  -  0  f 

The  r  1  diagrams  jre  ot  inlerest  from  several 
points  it  >iew  not  the  icasl  d  which  s  ihe 
demonstration  ot  the  imaging  soot  scattering  tech 
nic|ue  as  a  useful  probe  ot  tlame  structure  The 
interpolation  of  ihese  data  is  tairlv  straightfor 
ward  although  some  aveals  mav  tie  in  order  Ihe 
tusi  ot  ihese  is  that  the  interpretation  of  the 
absolute  intensity  ot  the  scattered  radiation  has  not 
been  established  As  was  mentioned  above  both 
•he  numher  lensitc  >t  panu  !es  mo  t tic  1  r  si/e 
lisiribution  determine  ihe  scattering  micnsiiv 
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Fig  4  Fig  ;  tor  ;  D  *  0  35 


Thus,  without  measuring  or  otherwise  evaluating 
one  of  these  parameters,  there  is  ambiguity  in  the 
interpretation  ot  the  intensity  An  unambiguous 
measurement  would  require  a  more  sophisticated 
technique  using  either  two-color  or  mulliangle 
scattering  at  each  pixel  location. 

The  most  interesting  and  instructive  aspect  of 
the  r-t  data  is  their  structure  As  will  be  argued 
below  this  structure  is  reflective  of  the  behavior  of 
the  (lame  sheet  surface  and  as  such  it  tells  much 
about  the  structure  of  the  flame  This  structure 
evolves  with  z/  D  and  is  quasi-penodic  in  time  at 
each  z  D  The  flame  structure  only  depends 
weakly  on  the  heat  release  rate  (fuel  flow  rate) 
over  the  range  studied,  even  though  the  overall 
length  of  the  flame,  Zt,  does  depend  on  the  fuel 
flow  rate  (see  Table  1) 

The  structure  of  the  soot  distribution  can  be  tied 
to  that  of  the  flame  sheet  by  several  arguments  In 
a  photographic  study  where  the  luminous  soot  and 


the  laser-illuminated  soot  were  observable  simul 
taneousiy  (data  not  shown;,  there  was  no  evidence 
of  laser-illuminated  soot  in  a  region  where  there 
was  no  luminous  soot  This  is  interpreted  to  mean 
that  there  is  no  appreciable  cold  (nonluminousi 
soot  in  the  flame.  This  suggests  that  the  soot  is 
located  primarily  close  to  the  hot  combusting 
flame  sheet 

Certainly  the  formation  of  soot  requires  the  high 
temperatures  found  near  the  diffusion  flame  If 
one  were  to  suppose  that  soot  were  to  survive  and 
move  away  from  the  flame  surface,  it  might  cool 
radiativeiy  and  show  up  as  nonluminous  soot— 
which  was  not  observed  Further,  and  perhaps 
more  significantly,  soot  surviving  away  from  the 
burning  regions  would  show  up  in  the  r-t  dia¬ 
grams,  coexisting  with  other  products  in  the  large 
heated  regions  of  the  flame  at  larger  z  D  stations 
This  is  also  not  observed,  nor  is  any  significant 
amount  of  soot  released  by  the  flame  as  a  whole 
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Fig  '  As  Fig  2  lor  ;  D  -  0  50.  except  here  ihe  vertical  lines  indicate  ihe  two  edges  ot  ihe 
burner  i  separated  hv  VI  ani  The  right  edge  <>t  (he  Figure  is  indistinct  due  lo  the  angular  culotl 
it  the  nterterence  tiller  The  angular  window  was  positioned  slightly  Dlfienter  to  ensure 
that  at  least  one  edge  of  the  flame  was  seen  clearly 


From  these  lines  of  evidence,  we  argue  that  for 
this  buoyant  natural  gas  flame,  the  soot  location 
approximately  marks  the  tlame  sheet  surface.  The 
correspondence  is  approximate  however  since  the 
soot  is  formed  on  the  fuel  side  of  the  tlame  sheet, 
before  the  fuel  encounters  oxvgen  and.  further,  the 
soot  hearing  region  is  probablv  larger  than  the 
reaction  rone  |4  10)  isee  Fig  12)  Thus,  the  soot 
distribution  is  probably  wider  than  the  tlame  sheet 
itself  in  a  direction  normal  to  the  flame  sheet 
surface  and  is  slightly  displaced  to  the  fuel  side  in 
the  same  direction  The  wider  soot  distribution 
will  not  misrepresent  the  structure  of  the  flame, 
however  it  the  s<x>t  bums  completely  before  it 
approaches  another  tlame  sheet  That  is.  the  s<x>t 
width  must  be  small  compared  to  the  Kolmogorov 
scale,  hut  not  necessarilv  small  compared  to  the 


tlame  sheet  thickness,  tor  the  stxit  distribution  to 
be  useful  as  a  marker  of  the  tlame  vheet  surtace 
The  difference  in  thickness  ol  the  two  surfaces 
must  nevertheless  be  kept  in  mind 

Thus,  where  there  is  soot,  there  is  a  tlame  sheet 
but  the  converse  is  not  necessarilv  true  It  is 
possible  that  in  an  appropriate  straining  field  ihe 
tlame  temperature  could  drop  below  that  required 
to  form  soot  but  not  so  low  as  to  extinguish  the 
tlame  { 1 1 1  Alternatively  the  effects  of  strain  mav 
he  to  diminish  the  time  available  tor  stxx  nuclea 
lion  and  growth  to  pr<x.eed  1 1 2 1  In  any  case  out 
data  do  not  allow  us  to  ascertain  whether  there  is 
combustion  without  stxx  formation  m  this  llamc 
hut  there  are  clear  indications  ot  breaks  in  the  soot 
distribution  tor  z  D  of  0  35  and  greater  The 
strain  rate  increases  with  the  velocity  as  c  P 
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increases  mi  this  may  possibly  be  viewed  as  a 
■.training  out  of  at  least  the  mhm  production 
Whether  the  entire  tlame  is  strained  out  at  the 
vame  places  is  not  addressed  bv  this  technique 
With  the  view  that  the  soot  structure  and  the 
tlame  sheet  structure  have  a  close  correspondence 
it  is  interesting  to  reconsider  Figs  Ml  Near  the 
burner  at  z  D  =  0  I  and  0  2  <  Figs  2  and  li  the 
structure  is  a  continuous  (involuted  surtace  The 
soot  thickness  which  is  presumably  an  upper 
bound  tor  the  tlame  sheet  thickness  as  discussed 
above  is  on  the  order  of  '  mm  Dm  estimate  was 
■  (burned  trom  the  minimum  widths  »t  the  mkk 
sheets  in  the  r  direction  in  the  r  i  images  At  about 
Z  D  =  0  15.  discontinuities  begin  to  appear  iFig 
4 1  and  the  structures  become  more  discontinuous 
with  increasing  ;  D.  although  the  spatial  resolu 
non  ot  our  measurements  mas  be  inadequate  at  or 
beyond  z  D  -  0  ’1  The  wmt  thickness  is  on  the 
order  ot  :  mm  at  :  D  -  0  15  The  spatial 


resolution  is  about  mm  so  hai  o  me 
thickness  approaches  mm  t  should  hfoms 
impossible  to  observe  all  'he  letaio  -t  tv  .... 
structure  This  is  consistent  auh'he  •bse'sco  a.  i 
ot  resolution  at  and  besond  ;  />  o'* 

One  inittails  surprising  'eaiurr  >t  he  .  ■ 
structure  is  that  the  .onioned  and  c.r  ar^e  /' 
broken  surface  does  mu  till  j  jriter  'racoon  t  hr 
space  Occupied  bs  ihe  danse  l*h, >ioki jphv  m,: 
rnosies  >t  ihe  uitunous  ..«.(  •  i.tuiion  • 

integrate  ail  uminous  mkh  ai.m^  'hr  me  i  -  ko  i 
lead  one  to  .onclude  *ronais  hai  he 
approaches  a  mote  >t  less  imtormn  mme 
filling  distribution  In  last  'he  .pposiir  .  see-' 
be  the  case  in  the  '  r  diagrams  *herr  be  ■•»  ' 
distribution  is  made  up  ilamelethke  'ragmrnr 
which  '^cupv  the  region  within  'he  >wlrr  n. >un,U 
ries  ot  the  tlame  thes  irr  tistnbme.i  ippi  . 
maiels  umlormls  >n  iseragr  hui  ml,,  dual.,  i" 
sharpls  localized 


LASER  S(X) r  SCATTERING  image 


1*1 


hr  mil  'hrti  *Oiu'  'I  'jnaturr  T  hr 

hIiui  ■<  ifi-Huir  n  m  'hr  Tiler  >l  'he  N.timo 
!"»•■<  rnjth  air  f  hr  Hn  soldi  sumhrr  'innn 
'  'T'  ihmii  ml'  Tear  hr  suinei  '*>  about  Nlli  seal 
! 1  ha\etj  'ft  rmlr  nlimatei  *i  he 

'iirr'  -ir  rt.<i  t>  hr  Mamr  iiametri  and ’he  h<e 
U"  •  rtrr'tjti  v  'tilt  T hr  '»rr»ll  '•nirlupe 
te*.  'I'tri  i  I'tpir  ;hi"  hape  ai  ihvi  ■  P  •  I 
*  irvi  h  i  ptjifin*  »  ir\  it  'hr  preitominani 
I  "•v»ur,t<  •  ’hr  :ui*>  ji  ji<r'  aiuri  T  •  P 
i'  ".*itr  i|  '  !  i  \«  •nr,nT*»iii  .<ii  *herti  inii'h 

T'l  '  <-*%  ifi.i.Tifii.  *  1 1 h i ft  hr  inline'  ’ 

•w  Til*  -hapr 

'Sr  mil'  'tint .  ’  hr  .*»x  ii«ir'hui«Ti  v  iris 

nmalcd  ■'  >  it<  •  in i  l  ippr r  unn  *hrtr 
i  'ixliii"S4  "»i  nirimtiirfH  •  >  pi<me«l  Thu 

•vtii">o*i  *»x  mr''»nlleiv  v  ■  Irlisril  n  hr 
sir  r'sillrsi  ’  hr  .*»x  ’return]  *'lhut  hr 
"hiihIjiT'  ’  ir  *h.  *ir  ''ante  I’  *  ai  as  mated 
'  >'•  lr>r"  11.11*  hr  ”.T'f  'nl4  den  *i  'hr 


>me  Jiiirihuium  n  'hr  ’  lata  set  S’i  ,-a«.h  'imr 
'nn*  iiaing  the  <>t*  threshold  *  alue  re  distinguish 
•he  rightmost  and  eftmuM  positions  ai  **h>ih  mux 
*41  pftirm  Then  'nr  r»».h  purl  'hr  'rastn's  .1 
iirte  'he  pnel  untamed  mi  *41  )is  tiled  hV  'hr 
’fa iiitm  >f  'ime  ihai  'hr  purl  *as  >niained 
between  the  rightmost  4ml  leftmost  uttnii  t  'he 
••»x  distribution  determined  >n  j  frame  ht  'rjirtr 
’vim  1'T  'he  alues  T  heal  'rlrjsr  ’jir  tno  •  P 
-sammed  herr  hr  . thIiIioiui  "in  nirmtutriv  . 
*  tratii  hr  'jittr  n  'hr  prn*xis,i  letned  '•»" 
nlrrmmrrH  i  teai  hr  rntr'  •!  hr  ’'ante  nv  . 
■hr  enter  >1  'hr  ’lame  i  n» m  *r  *  th">  hr 
'■lifn  'I  'hr  'lame  iivI  hr  Maine  -  .  xilimiai.* 
infirm  seat  'hr  rmrt  l-waid  'hr  -dir'  s.-«, 
Tft  he  Marne  -dgr  'tmth  it  and  'ut  *hrthr' 
hrrr  r  1 11I  InliiH'  put  f  i  t  h  X  .  and  'hr  xhI tit*  tij 
«*x  nirrmnierv  .  '  -uhsiantiaili  *>eaie'  ’han  hr 
simple  -i»x  hir'titiiirnv  *  li»r  hr  lull  *. 

'In'  ■<  hr  putlm*  ’rhii'.  t 


i 


0 


I 


|5 


10  2  0  3  0 

T  Mill) 

Fig  1 1  The  intensity  integrated  along  the  position  axis  is  plotted  as  a  function  of  lime  for  ;  D 
=  0.50.  The  full  3  5s  time  history  is  plotted  (with  every  two  frames  averaged),  showing 
approximately  seven  puffing  cycles 


Reaction  -  Soot  Laden 
Zone  Region 


Fig  12  The  left  sketch  shows  schematically 
tlamelet  The  right  sketch  indicates  the  regions 
a  tlamelet  is  swept  past  the  laser  beam 


Sheet 

location  of  soot  within  a  laminar  diffusion 
ere  scattered  laser  radiation  originates  when 


The  conditional  soot  intermittency  for  the  Z-  D 
=  0.2  case  is  plotted  in  Fig.  13  (upper  curve)  and 
is  nearly  constant  across  the  flame  There  is  an 
indication  of  a  dtp  near  the  center  i pixel  1 1 7)  as 
well  as  some  residual  structure  further  out  along 
the  radius,  but  the  curve  lies  predominantly 
between  45  and  55%  For  c  D  =  0  5  (Fig  14. 
upper  curve i.  the  conditional  soot  intermittency 
lies  near  40%  except  near  the  center  ipixel  260). 
where  the  peaks  to  about  70%  Figure  1  3  is  typical 
of  all  stations  less  than  ;  D  =  0  35,  while  Fig  14 
is  typical  for  axial  stations  z  D  =  0  35  and 
greater  The  peak  in  Fig  14  indicates  that  there  is 
some  difference  in  structure  between  the  center 
and  edges  of  the  soot  distribution  for  larger  z  D 


values,  which  we  believe  may  be  related  to  the 
connecting  strands  of  soot  between  putts  \e\er 
theless.  the  failing  off  of  the  simple  soot  intermit 
tency  toward  the  edges  is  not  present  in  the 
conditional  soot  intermittency  tor  any  of  ,iur  djta 
The  putf  shape,  then,  primarily  determines  be 
simple  soot  intermittency  curve  and  the  sond 
tional  soot  intermittency  is  generally  unitorm 
within  the  flame  edges  and  has  a  ialue  bear  -It 1  ' 
(neglecting  dtps  or  peaks  on  the  centerline •  The 
flame  sheets  mas  or  mas  not  be  js  discontinuous 
as  the  soot  sheets,  but  one  might  expect  that  be 
putf  shape  modulates  the  traction  of  time  'hat 
tlame  sheets  arc  encountered  also 

This  studs  raises  several  auestiv'ns  In  K 
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particular  tlames  studied.  the  soot  was  found 
primarily  in  thin  sheets.  The  entire  flame  was  in 
the  buoyancy  -dominated  range  and  was  a  inonpre- 
mixed)  diffusion  flume  If  one  is  led  to  surmise 
that  the  combustion  is  occurring  on  flame  sheets  in 
this  tlame.  a  supposition  that  perhaps  must  be 
supported  by  further  experimentation,  the  question 
arises  as  to  how  general  such  a  result  might  be  Do 
momentum-dominated  diffusion  flames  also  form 
sixn  sheets  and  does  the  combustion  -xcur  in  tlame 
sheets  in  these  tlames'  Greater  differences  might 
be  expected  with  premixed  combustion  since  a 
fuel  oxidizer  interface  does  not  exist  and  s<xn 
production  is  usually  quite  low 

To  answer  these  questions,  a  means  of  measur¬ 
ing  the  combustion  reaction  interlace  itself  must  be 
used  Two  dimensional  data  I r-t  or  r-o  must  be 
acquired  to  elucidate  the  structure  and  connectivity 
ot  the  reaction  zones  It  ma\  be  necessarv  to  work, 
n  a  nonsixnintf  tlame  ’such  as  hvdrogen  tlames i  in 
measuring  tlame  sheets  Jirectlv  to  avoid  interfer¬ 
ence  with  the  soot  particles  or  their  lurmnositv 

lONCl  l  SIGNS 

It  ta'  vn  shown  now  -he  scattering  ot  .user  ,ight 
•'v  s.H>t  t'urticies  uav  “x-  ;o  Jetermine  s.cm 

tistnbutions  n  jittusion  flames  'kith  same  ^ave 
its  be  st n it  ocatec  directed  adiacent  to  a  tlame 
-beef  .an  be  ’hough!  >t  as  marking  tlame  sheets. 
,i.tn.>u>th  ase'  mjv  ;x.st  where  a  flame  sheet 
s. .  jr-  p  ’be  absence  >f  adiacent  s.«ot  >r  n  other 
or  -'ames  >r  a  tb  '(her  fuels  mav  exist 

be  absence  ■'  i  •’urneiet 
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burner  is  also  consistent  with  an  increase  in  strain 
rate 

The  puffing  character  of  these  flames  was 
readily  apparent  in  the  data  and  was  found  to  be 
quasi-periodtc  with  a  frequency  of  about  0.5 
eg  D. 

Plots  of  soot  iniermittencv  indicate  that  the  soot 
is  present  most  frequently  near  the  burner  axis  and 
drops  to  zero  near  the  edges.  By  determining  the 
instantaneous  tlame  boundary  .  the  conditional  soot 
intermittency  can  be  analy  zed.  Plots  of  the  condi¬ 
tional  stxvt  intermittency  indicate  that  the  soot 
sheets  ix-oupy  about  40-60  “T  ot  the  instantaneous 
volume  within  the  outer  boundaries  of  the  tlame 
uniformly  jeross  the  flame  These  data  support  a 
mixiel  of  the  tlame  made  up  of  tlame  sheets 
distributed  approximately  uniformly  within  the 
confines  of  the  quasi-pertodic  pufflike  structures. 
This  mode!  indicates  that  the  surface  of  the 
stoichiometric  mixture  ratio  also  winds  throughout 
the  putts,  even  crossing  the  centerline 
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ABSTRACT 

Data  from  several  recent  experiments  on  mixing  and  chemical 
reactions  in  turbulent  shear  layers  and  Jets  is  discussed  in  some  detail 
and  used  to  formulate  a  picture  of  the  path  from  the  freestream  to  the 
molecular ly  mixed  state.  A  model  is  proposed  which  incorporates  the 
essential  steps  in  this  path  and  which  appears  to  provide  a  framework 
for  understanding  the  major  effects  of  Reynolds,  Schmidt,  and  Damktfhler 
numbers  on  the  chemical  reaction.  A  simplified  version  of  the  model 
reproduces  the  observed  Reynolds  number  dependence  of  nitric  oxide 
production  in  turbulent  fuel  jets. 


INTRODUCTION 


It  is  a  tenet  of  classical  turbulence  theory  that  in  the  limit  of 
Reynolds  number  approaching  infinity,  dissipation  is  independent  of 
Reynolds  number  and  takes  place  at  the  Kolmogorov  scale.  The  same 

assumptions  concerning  scalar  mixing  are  the  basis  for  a  model  of 
molecular  mixing  in  shear  layers  and  jets  that  has  been  inder 
development  for  3ome  time,  Broadwell  and  3reidenthal  (1982),  Broadwe.- 
and  Mungal  (1986,1987).  The  approach  was  suggested  by  the  experiments 
that  revealed  large  scale  organized  motions  in  shear  layers  and  jets, 
primarily  those  of  Brown  and  Roshko  (1974),  Konrad  (1976),  and 

Breidenthal  (1982)  in  shear  layers  and  of  Dimotakis,  Papantomou  and 
Miake-Lye  (1983),  in  jets.  In  the  shear  layer,  these  experiments  ana 

the  computations  of  Corcos  and  Sherman  (1976)  showed  that  instabilities 

lead  to  concentrated  regions  of  vorticity  and  that  free  stream  fluids 
enter  the  layer  in  streams  with  dimensions  of  order  of  the  layer 
thickness,  <5  .  It  is  a  basic  postulate  in  the  model  that  in  the  limit 
Re  -*  «  ,  no  mixing  takes  place  until  the  3cale  of  these  entering  streams 
is  reduced,  by  inviscid  motions,  to  the  Kolmogorov  scale,  \K  -  5/Re^/U  , 
where  Re  is  the  Reynolds  number.  When  the  scale  of  the  concentration 
fluctuations  in  the  streams  reaches  A  ,  it  is  shown  by  arguments  given 
later  that  the  time  required  for  diffusion  to  "homogenize"  the  mixture 
is  negligible  compared  to  the  time  to  reach  \K  from  5,  i.e.  compared  to 
<5/U  ,  where  U  is  a  characteristic  large  scale  velocity.  The  quantity, 
or  volume  fraction,  of  molecular ly  mixed  fluid  so  formed  is,  therefore, 
independent  of  both  the  Reynolds  number,  and  the  Schmidt  number,  Sc.  In 
this  limit  then,  the  rate  of  molecular  mixing  is  a  constant  as  is  the 
rate  of  entrainment  into  the  layer. 

When  the  Reynolds  number  is  finite,  but  such  that  (Re  Sc)  1/2  »,  , 
diffusion  layers  form  at  the  boundaries  of  the  entering  streams  and, 
because  their  thickness,  XT  ,  scales  with  the  large  scale  variables,  U 
and  <5  it  is  given  by 


1.  -  6/ (Sc  Re)1/2  -  6/Pe1/2 


wnere  ?e  is  a  3eciet  lumber.  Since  me  surface  irea  per  m.p  .  me, 

of  tnese  diffusion  sheets  scales  only  witn  o'  ,  their  fractiona.  /c . .m« 
3  •  at)  is  proportional  to  ?e*  2  . 

In  this  case,  also,  wnen  all  scales  in  tne  entering  stream  ire 
reduced  to  the  Kolmogorov  scale  the  mixture  is  Homogenized  is  oef:r». 
In  this  Lagrangian  description  of  tne  path  from  the  free  stream  •  :  me 
molecularly  mixed  state,  the  same  quantity  of  fluid  .3  .nvolved  ->r,ate /-m 
the  Reynolds  number.  Therefore,  the  fractional  volume  of  mixed 
produced  by  diffusion  at  the  Kolmogorov  scale  Is  independent  of  me 
Reynolds  number.  The  diffusion  at  the  scale  A,  causes  mixing  that  .a 
"early"  or  upstream,  in  the  tulerian  viewpoint. 

The  above  described  picture  of  shear  layer  mixing  is,  by  the  same 
argument,  applicable  to  Jets,  in  which  case  the  entering  stream  of 
reservoir  fluid  mixes  at  a  constant  rate  with  jet  fluid.  Differences 
between  the  two  flows  emerge  only  when  chemical  reactions  ire 
considered. 

At  any  axial  station  in  either  shear  layers  or  jets,  then,  the 
molecularly  mixed  volume  fraction,  Va,  can  be  written, 

Re  »  1 

vm  -  A  ♦  B/(Sc  Re)1  /2  (ScRe)1/2»i  (’) 

Re  »  (In  Sc)2 

in  which  A  is  the  average  volume  fraction  generated  when  the  entrained 
fluid  reaches  the  scale  and  B  is  the  average  non-dimensional  surface 
area  per  unit  volume,  i.e.,  S  measured  in  terms  of  5.  The  notation 
Sc  Re  is  used  instead  of  Pe  because  of  the  need  to  state  the  conditions 
of  applicability  of  the  model.  The  origin  of  the  restriction,  Re  » 
(In Sc)2  is  discussed  later. 

In  the  above  referenced  discussions  of  the  model,  the  diffusion 
layers  were  called  strained  laminar  flames.  Connection  is  made  with  the 
preceding  discussion  by  noting  that  if  these  layers  are  in  equilibrium 


wi  tn  '  ,;ci.  .  irge  r.'i.".  *  -  •  ,  mem  m.mr.ess 

propor  t .  ma  1  to  2  t  "  and  me  expression  above  -esui'.s.  j«e  .ar '•.•*'■. 
Fende..,  and  Marti*  '  for  a  Jisooaa.tn  of  mese  f.ames.  3  .3  ■  -e 

JiffiSion  :cef  f  id  i*ant .  When  ic  *  ♦ ,  is  me  ’iy.tr  soa.*?.  i  '  i  •• 

brought  to  our  attention  oy  -(.  4.  ...epmann  wno  points  tut  tnat  *”.e  . 

are,  in  the  aathemat  I  cal  sense,  Interna,  ooundary  layers. 

The  Jiscussion  so  far  na3  oeen  toncerned  tn.y  with  •no.ecu.ar  -,.<m.g 
with  no  restriotions  tn  the  nature  of  the  tnemioa.  "eact.tns. 
particular,  reaction  rates  of  any  magnitude  can  5e  considered.  Th.s 
flexibility,  also  a  character istic  of  the  coherent  flame  model,  ‘Martin 
and  Broadwell  ( 1 977 ) ,  may  be  useful  in  interpreting  the  -scent ly 
measured  finite  kinetic  effects  in  hydrogen  flames,  Dibble  and  ^agre 
(1987),  and  in  studying  lift-off  and  Plow-out  of  fuel  Jets. 


With  regard  to  the  Taylor  diffusion  layers,  notice  that  in  the  jet 
they  form  between  newly  entrained  reservoir  fluid  and  jet  fluid,  i.e. 
fluid  that  is,  or  was,  "turbulent".  For  them  to  persist  it  appears  to 
be  necessary  for  the  small  scale  turbulent  motions  in  the  jet  fluid  to 
have  been  destroyed,  by  momentum  diffusion,  when  the  concentration 
variations  were  destroyed.  In  the  shear  layer,  this  picture  is 
consistent  with  the  observations  of  Coles  (1983)  and  Hussain  (1983), 
that  turbulence  is  generated  in  the  braids  and  dissipated  in  the  cores 
of  the  vortices.  Furthermore,  if  the  small  scale  motion  is  dissipated, 
the  vortices  would  be  left  in  solid-body-like  rotation,  evidence  for 
which  is  reported  by  Wygnanski  and  Fiedler  (1979).  In  both  shear  layers 
and  jets,  the  motions  are  unsteady  on  the  largest  spatial  and  temporal 
scales  in  the  flow.  These  unsteady  motions,  as  illustrated  by  the  flame 
length  fluctuations  observed  by  Dahm  and  Dimotakis  (1985),  for  instance, 
make  possible,  if  not  plausible,  the  idea  that  the  small  scale  motions 
and  concentration  fluctuations  are  successively  generated  and  destroyed 
in  the  Jet  fluid  as  it  moves  along  the  axis. 


While,  as  Just  stated,  large  scale  unsteadiness  is  an  Inherent 
feature  of  the  flows  being  considered,  the  model  attempts  to  deal  only 
with  their  averaged  consequences.  In  addition,  in  most  of  the 
discussion  that  follows,  the  density  will  be  taken  to  be  constant.  The 
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THE  HIGH  RffTMOLOS  NUWER  LIMIT 


In  the  above  cited  references  to  earlier  worn  on  the  model,  an 
upper  limit  for  the  time  required  for  diffusion  to  homogenize  the  fluids 
at  the  Kolmogorov  scale  was  estimated  from  the  time  to  diffuse  across 
that  scale,  l.e.,  from  1*^/0  -  (R/u)  Sc/Re1 /2  where  R  is  the  jet  radius 
and  u  the  axial  velocity.  Batchelor  (1959)  shows  that  local  straining 
reduces  X<  to  R/(Re^/USc1  /2)  in  a  time  (R/u)  (l/Re1/2)  In  Sc  ,  a  result 
brought  to  the  author's  attention  by  P.  E.  Dimotakls.  Since  the  time 
for  diffusion  across  the  smaller  scale  is  only  (R/u)/Re1/2  ,  the  larger 
time  (i/Re1/2)ln  Sc  is  controlling  and  is  the  restriction  given 
relative  to  Equation  (1). 

When  the  conditions.  Re  »  1,  (Sc  Re) 1/2  »  1,  and  Re  »  (In  Sc)2  are 
satisfied,  it  is  a  simplification  to  consider  the  scales  to  be  divisible 
into  the  three  separate  categories: 

R,  R/Re1 /2  and  R/Re3/4  , 

and  to  assume  instantaneous  mixing  at  the  scale  XK. 

Consider  first  the  limiting  case  Re  ♦  »  for  any  fixed  Schmidt 
number,  in  which  case  mixing  takes  place  only  at  \K.  (Note  that  the 
requirement  Re  >  (In  Sc)2  precludes  the  case  Sc  ♦  •  for  fixed  Reynolds 
number.  It  seems  clear,  intuitively,  that  the  mixed  volume  fraction 
goes  to  zero  in  this  limit.) 

A  critical  check  on  the  model  in  the  high  Reynolds  number  limit  is 
whether  or  not  the  mixing  rate  is  Independent  of  Re  and  Sc.  Very  little 
evidence  is  available  concerning  this  point,  but  the  experiments  of 


3re i ientha  1 


’  98’  and  Koocnesf anan i  and  UimotaKis  '966 


:n  tne  same 


facility  and  over  a  limited  Reynolds  number  range,  "O4  <  Re  <  ’  D3  ;  snow 
that  for  3c  '  oOO,  tnere  is  no  dependence  on  Re.  These  results  suggest 
tnat  3c  Re  '  ’ 3  .s  sufficient  for  tne  limiting  condition  to  app.y.  In 
gases,  Sc  '  ’  ,  at  Re  "  '0^  ,  nowever,  Mungal,  Hermanson,  and  Dimotarfis 
'.'905;  find  Reynolds  number  effects;  '  Q3  is  apparently  not  nign  enoug.n. 

Indirect  support  for  tne  proposed  mixing  mechanisms  m  tne  nign 
Reynolds  number  limit  is  in  tne  composition  and  spatial  distribution  of 
tne  molecularly  mixed  fluid.  If  there  is  a  delay,  after  entrainment, 
during  which  the  fluids  are  distributed  by  large  scale  motion  throughout 
the  shear  layer  or  the  Jet,  the  mixed  fluid  composition  should  be 
independent  of  the  lateral  coordinate  in  shear  layers  and  of  the  radius 
in  axi-symoetric  Jets. 


Konrad  (1976)  noted  the  approximate  flatness,  in  his  gas  shear 
layer,  of  the  mixed  fluid  composition  and,  in  fact,  this  observation  was 
one  of  the  motivations  of  the  present  model.  The  most  remarkable 
results,  however,  are  those  obtained  in  shear  layers  in  water  (Sc  -  600) 
and  Re  -  2  •  i  04  by  Koochesfahani  and  Dlootakis  (1986). 

Figure  1  shows  their  mean  mixed  fluid  composition  deduced  from  chemical 


reaction  experiments  that  eliminated  probe  resolution  errors.  (They 
show,  incidentally,  the  rigorous  requirements  for  accurate  measurements 
of  this  quantity  by  conventional  Instrumentation).  We  see  that  there  is 
no  indication  whatever  of  a  lateral  gradient.  It  is  a  convenient  and 
useful  additional  approximation  for  modeling,  used  the  following,  to 
assume  that  the  mixed  composition  has  a  single  value,  but  what  is  to  be 
noted  here  is  lateral  distribution.  It  is  difficult  to  escape  the 
conclusion  that  the  fluids  are  distributed  throughout  the  layer  before 
they  mix. 

The  picture  is  not  as  simple  for  Jets.  Dahm  and  Dimotakis  (1985), 
again  in  experiments  in  water  in  the  Re  range  103  -  101*,  conclude  from 
laser  sheet  photographs  and  composition  profiles  that  there  are  large 
region  of  nearly  uniform  composition  but  that  they  are  so  arranged  that 
radial  scans  sometimes  Intersected  more  than  one.  Figure  2  contains 
samples  of  this  measurement.  More  data  is  needed  but  the  few  results 


Figure  1.  Dependence  on  lateral  coordinate,  y/6,  of  avenge 
concentration  of  mixed  fluid.  CA  denotes  high 
speed  fluid.  Koochesf ahani  and  Dimotakis  (1986). 

are  In  accord  with  the  notion  that  large  scale  motions  distribute 
reservoir  fluid  throughout  the  jet  before  it  mixes  molecularly.  Such 
motions  are  dramatically  illustrated  in  high  speed  laser  sheet  images  of 
Long  (1987).  Likewise,  temperature-time  traces  in  a  plane  Jet,  Antonia 
et  al  (1986),  have  been  shown  by  them  to  be  evidence  of  large  scale 
motions. 


JET  MODEL  OF  THE  LIMIT  Sc  Re  ♦  •  ,  Re  »  1  ,  Re  »  (In  Sc)2 

To  put  the  above  description  in  approximate  quantitative  form,  we 
consider  only  the  far  field  where  the  jet  spreads  linearly  and 
similarity  in  the  velocity  distribution,  u  (x/dQ,  n)  has  been  attained. 


Here,  x  :a  the  jx:a.  cistance,  -re  et  "c::.e  ccimerer,  3,  - 

Jet  radius,  and  n  •  r  R  .  L-iwewiae  the  mcleca-ar  ly  n.  xed  f-a ct. 
ia  taxen  to  depend  only  on  n  .  Recall  tnat  tne  juat.f;.  cat .  :r.  : 


independence  of  a  from  x  la  tnat  tne  entrained  flute  -eaones  tne 
Kolmogorov  acale  and  mixes  after  it  moves  downstream  a  distance 
scales  only  with  tne  -adius  at  wnion  it  was  entrained. 

The  connection  between  a  and  tne  rate,  /*,  at  wn.cr.  "“servo."  :  ' .  _ . 
"enters"  tne  moiecularly  -nixed  flax  is 


d  p 

vh  »  2ir  —  R  i  a(  n)  a  v  X/  dQ  ,  n ;  n  dn 
dx 


R  -  Ci1/r2x  and  u  -  (uq^/x)  g(n) 


Then, 


where  E 


m  2-n  C,  u0  d0  E 


a(  n)  g(  n)  n  dn  . 


The  Independence  of  a  from  x  implies,  therefore,  that  the  molecuiarly 
mixed  fluid  flux  rises  linearly  with  x  as  does  the  total  Jet  flux.  As 
ta  well  known,  the  latter  result  follows  simply  from  the  conservation 
equations;  the  numerical  rate  has  been  measured  by  Ricou  and  Spalding 
(1961 ). 


In  the  shear  layer,  as  was  noted  above,  the  mean  concentration  In 
this  limit  is  independent  of  the  lateral  coordinate.  If,  in  addition, 

the  concentration  is  taken  to  have  a  single  value,  then  the  integral 

► 

o(n)  dn  ,  the  average  mixed  volume  fraction,  can  be  determined  from 
chemical  reaction  experiments  to  be  0.28.  (See  Broadwell  and  Mungal, 
1986.)  Dahm  (1985)  computed  the  mean  mixed  concentration  from  profiles 
like  those  in  Figure  1,  finding,  of  course,  a  radial  variation  in  this 
quantity.  Figure  3*  Such  variation,  as  well  as  probe  resolution 
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A  *  a  i;  ii  i  >  0.5. 

This  value,  twice  tnat  in  tne  shear  layer,  is  in  approximate  accord 
with  that  found  by  Waddell  in  his  use  of  acid-base  reactions  in  water  to 
study  mixing  in  Jets,  Hottel  (1953)-  (This  early  work  provided  tne 
starting  point  for  the  liquid  shear  layer  and  jet  experiments  discussed 
above. ) 

With  the  assumption  that  tne  mixed  fluid  composition  is  independent 
of  n  and  has  a  single  value,  the  fuel  conservation  equation  becomes 

2m—  (Ce)  R1  I  a(  n)  u(x,  n)  ndn  «-2mk(Cp)  (Co)  R2ia(n)-:dn 
dx  r  h  J  r  h  u  h 

The  oxidizer,  carried  in  the  reservoir  fluid,  is  governed  by 


2”S(C0)h«’ 


a(n)u(x,n)  ndn  •  -2m  <  (Cpl 


-O', 


where  (Cp)^  and  (C(j)h  denote  the  fuel  and  oxidize'-  concert -a:  . 
mixed  fluid,  (Cq).  is  the  reservoir  oxidizer  concer.t-at  . -r .  r 
reaction  rate  coefficient  for  the  reaction: 
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binary  reaction  with  a  constant  rate  coefficient.  As  will  be  seen,  this 
simpif icatlon  can  be  dropped  when  numerical  computations  and  complex 
reactions  are  considered. 


Introducing  the  variables  C*0  .  (Co)h/(CQ).  ,  C*F  -  (CF)h/ (CF) i  , 
and  z  ■  x/d0  ,  where  (Cp) ^  in  the  initial  fuel  concentration,  and  using 
Equations  (2)-(5),  we  can  write  Equations  (7)  and  (8)  in  the  form 


—  (zC»F)  -  -Da  z*  C*F  C*0 
dz 


—  (2C*0) 
dz 


-Da  $  z*  C*p  C*q  *  t 


in  which  4  is  the  equivalence  ratio,  (Cp)i/(CQ)-,  and  Da  the  Damkohler 
number  defined  by 

A( Cq ) m  k  dQ 


and  A  -  o ( n)  n  dn  . 


When  there  is  no  reaction.  Da  »  0,  Equation  (9)  yields 


C*F  -  c/z  -  3.2/z  01) 

where  the  constant  is  taken  from  Equation  (6). 

Since  (C*F  +  c*h)  "  1  for  Da  -  0  , 

C*F  -  1  -  3.2/z 

Multiplication  of  Equation  (9)  by  $  and  subtraction  from  Equation 
(10)  yields 


—  [z  (C*0  -  *  C»F]  -  1  . 
dz 


from  which  we  find. 


<C*o  “  ♦  C  p) 


Z-3.2  (♦♦! ) 


in  which  the  constant  comes  from  Equation  (11).  For  Da 
the  flame;  therefore 
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and  the  flame  ends,  C*g.  .  o,  at 

ze  -  3.2  U-1  )  (12) 

Since  Equation  (12)  is  valid  for  any  form  of  the  reaction  term,  as 
long  as  Da  •*  •,  it  may  be  compared  with  the  "flame  length"  measurements 
of  Waddell  and  of  Dahm  et  al  (1  984)  for  reactions  in  water  where  Ra  * 
103  -  104  and  Sc  Re  "  106  -  107. 

The  experimental  results  are  shown  in  Figures  4  and  5  ,  the  first 
establishing  the  Independence  of  ze  from  Reynolds  number  for  Re  >  3000  , 
and  the  second  the  linear  dependence  of  ze  on  <p.  Since  in  the 
experiments  zg  was  taken  to  be  the  axial  location  of  the  last  indication 
of  any  unreacted  Jet  fluid  and  since  the  flame  length  fluctuations  are 
large,  it  is  understandable  that  the  measured  lengths  are  much  larger 
than  that  given  by  Equation  (12).  The  important  point,  however,  is  that 
the  linear  dependence  of  ze  on  <p  in  the  experiment  is  consistent  with 
the  Equation  (12),  and  hence  consistent  with  the  assumption  that  a  is 
independent  of  axial  distance. 


FINITE  REYNOLD  NUWER  JET 


At  Reynolds  numbers  for  which  the  Jet  is  turbulent  but  low  enough 
for  the  Taylor  layers  or  strained  flames  to  make  significant 
contributions  to  the  molecular  mixing,  the  procedure  is  as  follows.  To 
simplify  the  notation,  the  Schmidt  number  will  be  taken  to  be  unity  in 
this  section. 


•rv 


Recognizing  that  the  strained  flames  form  between  the  newly 
entrained  reservoir  fluid  and  the  continuously  reacted  jet  fluid,  we  may 
sketch  the  process  as  shown  in  Figure  6.  The  equality  of 


Figure  6.  Schematic  diagram  for  reactions  in  the  jet  model. 


the  net  volume  flux  to  the  flame  sheets,  v^,  from  the  homogeneous 
mixture  and  the  free  stream  comes  from  the  assumption  of  constant 
molecular  transport  properties,  an  assumption  leading  to  sheets 
constituted  of  equal  amounts  of  fluid  from  the  two  streams.  Motice  that 
the  net  flux  to  the  homogeneous  mixture,  vhf  i3  unchanged;  the  presence 
of  the  flame  sheets  only  changes  the  constituents  in  the  entering 
stream.  The  fluid  in  the  flame  sheets  at  any  axial  station  becomes  part 
of  the  homogeneous  fluid  further  downstream. 

The  net  volume  addition  to  the  flame  sheets,  vf,  can  be  written 


2x 


5"  f 


B(n) 

R 


b  R 


Re 


1/2 


u  n  d  n 


(13) 


in  which  0(n)/R  is  the  flame  sheet  area  per  unit  volume  and  the  sheet 
thickness  is  taken  to  be  that  of  a  strained  flame  in  equilibrium  with 
the  large  scale  strain  u/R  ,  i.e. 


u(  i.t  i.M.t  i.i  i>  < 


XT  -  b  (D/c)1/2  ■  b  (DR/u)1/2  -  b  R/(Re)1/2 


where  D  Ls  a  constant  of  order  one  (See  Carrier,  Fendeil,  and  Marble, 


1975.)  With  Equations  (3),  Equation  0  3)  becomes 


'f  -  2 it  Cl  u 0  d,  F/Re1 


where 


E  *  &  3( n)  g(n)  n  d  n  . 


As  before 


vh  -  2ir  Cj  u0  d0  E 


E  *  |a(n)  g(n)  n  d  n 


With  the  help  of  Equations  (3)  and  (14),  fuel  and  oxidizer 
conservation  equations  similar  to  Equations  (9)  and  (10)  can  be  written 
for  both  the  flame  sheets  and  the  homogeneous  mixture.  Now,  however, 


the  flame  sheet  chemical  reactions  must  be  treated  numerically  even  for 


the  simple  reaction  considered  in  the  preceding  section.  The  proposed 
scheme  for  solution,  currently  being  carried  out  for  a  simple  reaction, 
is  to  consider  at  each  integration  step,  the  flame  sheet  properties  to 
be  those  of  a  strained  flame  in  equilibrium  with  the  local  strain,  u/R, 
and  to  use  either  numerical  solutions,  such  as  those  of  Dixon-Lewis  et 
al  (1984),  or  approximations  to  them.  One  side  of  the  flame  is  pure 
reservoir  fluid  and  the  other  the  fluid  of  changing  composition  from  the 
homogeneous  mixture. 


In  the  absence  of  such  solutions,  however,  several  useful 


conclusions  can  be  drawn  from  the  nature  of  the  model.  First,  since  the 


flame  ends  when  all  the  fuel  is  consumed  in  both  regions,  and  since  the 


addition  rate  to  the  homogeneous  mixture  is  unaffected  by  the  flame 
sheets,  the  flame  length  is  independent  of  Reynolds  number  in  this  case 
also.  The  same  argument  shows  the  independence  from  Schmidt  number. 
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Experimental  information  on  this  point  is  discussed  later. 

Julian  Tishkoff  called  this  paradox  to  our  attention,  i.e.,  the 
apparent  contradiction  between  thf  \  "luence  of  Reynolda  and  Schmidt 
numbers  on  molecular  mixing  and  the  lependence  of  flame  length  from 
these  parameters.  The  arguments  just  given  appear  to  resolve  the 
paradox:  the  upstream  molecular  mixing  in  the  strained  flames 
influences  the  mixing  rate  but  does  not  reduce  the  distance  required  for 
every  element  of  Jet  fluid  to  be  mixed  with  $  parts  of  reservoir  fluid. 

Next,  the  model  provides  a  basis  for  understanding  the  dependence 
of  the  nitric  oxide  production  rate  on  Reynolds  number.  The  experiments 
of  Bilger  and  Beck  0  974)  on  hydrogen-air  flames  and  of  Peters  and 
Donnerhack  (1901)  on  methane-air  showed  that  when  the  effect  of 
residence  time  was  removed,  there  remained  an  effect  of  Reynolds  number; 
the  amount  produced  appeared  to  decline  approximately  with  Re'/2. 
Since  NO  is  produced  by  a  slow  chemical  reaction  when  the  temperature  is 
high,  in  the  model  it  appears  in  the  flame  sheets  only  at  small  x/d0  , 
where,  in  the  absence  of  radiation,  the  adiabatic  flame  temperature  is 
reached.  The  forward  NO  reactions  stop  and  NO  may  be  destroyed  when  the 
flame  sheet  material  returns  to  the  homogeneous  mixture.  Near  the  end 
of  the  flame  NO  is  produced  in  both  regions.  As  Equation  '14) 
illustrates  explicitly,  the  flame  sheet  NO  production  also  depends  on 
1 /Re1 /2  and  the  model,  therefore,  has  a  mechanism  for  explaining  the 
experimental  results. 

To  examine  these  ideas  in  more  detail,  an  approximate  calculation 
was  carried  out  for  both  fuels  with  a  full  set  of  kinetic  equations.  A 
representative  list  of  the  reactions  and  species  can  be  found  in  Heap, 
et  al  0  976),  and  the  actual  updated  set  used  in  the  calculations  in 
Bowman  and  Corley  (in  press). 

For  the  model  computation,  Tyson,  Kau  and  Broadwell  0  982),  the 
strained  flame  "reactor"  was  considered  to  contain  a  homogeneous  mixture 
like  that  generated  at  the  Kolmogorov  scale,  but  vf  from  the  homogeneous 
mixture  was  varied  to  maintain  a  stoichiometric  mixture  in  the  flames 
and  thus  to  yield  the  adiabatic  flame  temperature  reached  in  strained 


flames.  The  constant  fixing  the  flow  rate,  vh>  was  inferred,  as 
described  above,  from  the  flame  length  data  of  Waddell  (seee  Hottei 
1953)  and  the  remaining  constant  chosen  to  match  the  model  to  the  data 
at  one  point.  Approximate  account  was  taken  of  radiation  in  both 

flames. 

The  model-experiment  comparison  for  the  hydrogen  flame  ;s  shown  in 
Figure  7.  Only  the  NO  emission  index  was  measured  for  the  methane  flame 
and  while  the  agreement,  with  the  same  constants,  in  that  case  ;s 
equally  good,  the  Reynolds  number  range  for  the  high  Froude  number  runs 
was  limited.  In  fact,  the  Reynolds  numbers  of  both  experiments  are  30 
low  that  the  good  quantitative  agreement  between  the  model  predictions 
and  the  experiments  is  likely  only  fortuitous.  The  qualitative  trend 
with  l/Re1^,  however,  is  clear  in  the  model  and,  having  been  observed 
in  different  experiments  in  two  gases,  seems  believable.  There  are 
plans  to  repeat  the  experiment  and  to  make  a  more  complete  model 
calculation. 

A  more  straight-forward  check  of  the  ability  of  the  model  to  deal 
with  molecular  effects  is  a  comparison  of  model  predictions  with 
experiment  of  the  differences  between  gas  and  water  reactions,  Mungal 
and  Dimotakis  (1984).  The  satisfactory  agreement  in  that  case  as  well 
not  only  gives  more  confidence  in  the  model  but  lends  credence  to  the 
idea  that  nitric  oxide  production  in  flames  is,  indeed,  Reynolds  number 
dependent . 


EFFECTS  OF  DENS ITT  CHANGE 


In  the  discussion  so  far  the  density  has  been  taken  to  be  constant. 
However,  Ricou  and  Spalding  (1961)  show  that  when  the  Jet  density,  oi( 
differs  from  that  In  the  reservoir,  p„f  the  ratio  of  the  total  Jet  flux 
to  the  Initial  value  Is 
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Figure  7.  Comparison  of  calculated  nitric  oxide  concentration  with 
measurements  for  hydrogen-air  flame.  Q .  Re  "  1.540; 
A,  8e  >  4,350;  O.  Re  *  12,300.  (Froude  number 
'■6-10  .  Data  from  Bilger  and  Beck  (1974).  From 
Tyson,  et  al.  (1982). 
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Since  this  relation  is  consistent  with  the  idea  that  the  entrainment  is 
controlled  by  the  jet  momentum,  it  may  be  expected  to  hold  when  the  jet 
is  reactive.  Avery  and  Faeth  (1974)  develop  a  theory  for  such  jets  and 
with  it  are  able  to  correlate  the  flame  lengths  of  a  remarkable  variety 
of  reactants  by  a  relationship  of  the  form  of  Equation  (15)  in  which 
is  replaced  by  d*-d0(p„/p  j ) 1 /2  .  Dahm  it  Dimotakis  (1985)  present  tr.eir 
data  together  with  that  collected  by  Avery  and  Faeth  in  this  form  is 
shown  in  Figure  8. 

The  results  in  this  figure  are  explained  by  the  present  model  if 
the  distance  after  entrainment  for  the  jet  and  reservoir  fluid  to  reac.n 
the  Kolmogorov  scale  continues  to  depend,  in  this  case  also,  only  on  the 
local  jet  radius.  The  arguments  leading  to  Equation  (11)  appear  to  oe 
valid  in  this  more  general  situation. 


DAMJ«3hL£R  number  effects 

Experiments  by  Mungal  and  Frieler  (1984)  on  the  hydrogen- fluorine 
reaction  is  a  shear  layer  cover  the  Damktfhler  number  the  range  from  zero 
to  mixing  limited  conditions.  In  this  flow  and  for  their  fixed  high 
equivalence  ratio,  an  integral  treatment  of  the  strained  flame  allows  a 
complete  analytical  solution  of  the  model  equations.  Figure  9  shows  a 
comparison  of  the  solution  with  the  experimental  results  (Broadwell, 
Mungal,  1986).  The  constants  A  and  B  in  the  model  equations  had  been 
determined  from  earlier  experiments  (Mungal  and  Dimotakis,  1984)  on  the 
fast  reaction  in  this  flow  and  from  the  experiments  in  water 
(Koochesfahani  &  Dimotakis,  1986).  In  the  figure  Da  is  defined  as 
k(Cjj)-  x/u  where  k(C^)«  is  the  reaction  rate  coefficient,  x  the  distance 
of  the  measuring  station  from  the  splitter  plate,  and  u  the  mean  shear 
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Figure  9.  Dependence  in  shear  layer  of  product  on  Damkohler  number, 

Da  *  x/u.  6p,  total  product  thickness;  (5p)f,  flame 

sheet  product;  (<5p)h>  homogeneous  mixture  product;  5,  shear 
layer  thickness.  The  symbols  are  experimental  data  from 
Mungal  and  Frieler  (1985)  for  two  different  velocities  at 
the  same  velocity  ratio. 
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layer  velocity.  The  ordinate  fip/5  is  a  non-dimensional  measure  of  tne 
amount  of  product  at  this  station.  The  separate  contributions  of  tne 
flame  sheets  (<Sp)f  and  the  homogeneous  zones  (50)^  are  shown,  with  the 
latter  seen  to  reach  its  limiting  value  more  quickly. 

In  the  jet,  the  varying  jet  fluid  composition  appears  to  preclude  a 
similar  integral  treatment  of  the  flame  sheets.  Hence,  only  the  nigr. 
Reynolds  number  situation,  described  by  Equations  (9)  and  (10),  nas  seen 
considered.  C.  F.  Frieler  (private  communication)  notes  that  an 
analytical  solution  of  these  equations  is  possible  but  that  it  is  of 
such  complexity  that  numerical  results,  which  he  has  also  kindly 
provided,  are  more  useful.  These  are  shown  in  Figure  10,  in  whicn  we 
see  that  the  limiting  concentration  distribution  is  attained  for  0a  >  ’ 
compared  to  Da  >  10  in  the  shear  layer.  The  decline  in  jet  velocity 
with  distance  must  be  responsible  for  this  difference  as  well  as  for  tne 
non-linear  influence  of  Da  on  the  jet  solution.  The  flame  length  for 
Da  -*  «  ,  i.e.  >  1  ,  is  32  in  agreement  with  Equation  (11). 


CONCLUDING  REMARKS 


As  has  been  pointed  out  in  the  foregoing,  an  implication  of  tnis 
analysis  is  that  extremely  high  values  of  the  parameter  Sc  Re  are 
required  to  eliminate  explicit  molecular  effects  on  molecular  mixing. 
For  a  Schmidt  number  of  unity,  it  appears  that  Reynolds  numbers  of  the 
order  of  10®  -  10?  are  necessary.  This  would  mean  that  most  laboratory 
combustion  experiments  are  significantly  influenced  by  "low  Reynolds 
number"  effects  and  that  analyses  neglecting  them  are  likely  to  be 
inadequate.  These  observations  may  also  have  some  bearing  on  the 
discussion  of  fluid  mechanics  at  high  Reynolds  number  as  noted  in  the 
Introduction. 
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With  regard  to  flame  sheets,  the  present  analysis  emphasizes  tne 
notion  that  the  fuel  side  consists  of  fuel  mixed  with  a  linearly  rising 
amount  of  products,  from  x/dQ  -  10,  and  that  the  flame  sheet  concept  is 
appropriate  for  any  DamkBhler  number  provided  the  Peclet  number  is 
sufficiently  high.  The  proposed  model  is  simple  enough  that  output  from 
fairly  elaborate  strained  flame  calculations  should  be  able  to  be 
incorporated  fairly  easily.  Since  there  is  no  limitation  on  the 
DamkOhler  number,  the  same  sort  of  calculations  may  help  clarify  the 
mechanisms  controlling  lift-off  and  blow-off. 
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